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construction  period.  Objectives  of  this  research  program  are  to  develop  tech¬ 
niques  to  minimize  anil  control  scour  during  nearshore  construction,  and  to 
predict  the  probable  magnitude  of  scour  that  may  result  as  a  function  of  cur¬ 
rents  and  wave  climate.  One  phase  of  the  research  program  is  development  of 
numerical  techniques  (incorporating  both  refraction  and  diffraction  effects 
near  the  structure)  for  computing  wave-induced  velocities,  tidal  currents,  and 
wave  heights  in  the  vicinity  of  structures,  and  applying  these  results  to  de¬ 
termine  sediment  transport  of  the  bottom  material  at  the  particular  site. 

The  present  state  of  nearshore  current  and  wave  theories  has  reached  the 
point  where  detailed  experimental  investigations  are  required  for  the  verifi¬ 
cation  of  analytical  developments  and  numerical  models.  To  provide  a  founda¬ 
tion  for  further  advancements,  a  simple  beach  profile  consisting  of  straight, 
uniform  contours  parallel  with  the  shoreline  was  experimentally  studied  by 
Hales  (1980).  A  shore-connected,  vertical,  thin,  impermeable  barrier  (break¬ 
water)  was  installed  perpendicular  to  the  shoreline  to  simulate  prototype 
jetties  and  breakwaters  commonly  occurring  along  many  coasts.  The  purpose  of 
the  present  study  is  to  extend  the  previous  work  of  Hales  (1980)  by  installing 
a  shore-connected,  vertical,  thin,  impermeable  breakwater  at  a  60-deg  angle  to 
the  shoreline  to  simulate  a  larger  range  of  prototype  jetties  and  breakwaters 
in  existence  at  the  present  time.  Experimental  measurements  of  refraction  and 
diffraction  downcoast  ol  this  oblique  structure  were  made  to  obtain  quantita¬ 
tive  knowledge  of  this  phenomenon  in  the  lee  of  the  .jetty  or  shore-connected 
breakwater.  These  data  weie  then  compared  with  the  uniformly  valid  asymptotic 
theory  of  l,iu,  I.ozano,  and  Pantazaras  (1979  )  for  the  same  arrangement. 

The  numerical  model  for  determining  wave  heights  downcoast  of  a  straight 
breakwater  at  an  angle  to  the  shoreline  under  combined  refraction  and  diffrac¬ 
tion,  based  on  the  uniformly  valid  asymptotic  theory,  was  obtained  by  contract 
with  Dr.  Philip  Liu,  Cornell  University.  Because  the  uniformly  valid  asymp¬ 
totic  theory  is  developed  from  the  small  amplitude  wave  assumptions,  the  ef¬ 
fect  of  varying  the  incident  wave  height  on  the  wave-height  amplification  fac¬ 
tor,  H/H  ,  was  investigated.  The  theory  and  experimental  data  were  found  to 
compare  favorably  tor  all  wave  heights  tested  in  the  shadow  zone,  but  diverged 
with  increase  in  incident  wave  height  farther  downcoast  out  of  the  shadow 
zone.  While  the  uniformly  valid  asymptotic  theory  is  far  superior  to  diffrac¬ 
tion  theory  alone  under  these  conditions,  additional  numerical  work  should 
incorporate  nonlinear  wave  theories  for  completeness. 
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Assistant  Chief  of  the  Hydraulics  Laboratory;  R.  A.  Sager,  Clue!  nt  tin  !  •.lu¬ 
rries  Division  and  IOMT  Program  Manager;  Dr.  R.  W.  U'halin.  former  Chief  I  'he 
Wave  Dynamics  Division;  Mr.  D.  D.  Davidson,  Chief  of  the  Wave  Reseurih  Hr  iin  li; 
and  Dr.  J.  R.  Houston,  Research  Engineer  and  Principal  Investigator  tor  the 
Erosion  Control  of  Scour  During  Construction  work  unit.  The  Wave  Dynamics  Di¬ 
vision  was  transferred  to  the  Coastal  Engineering  Research  Center  (CERC)  of 
WES  on  1  July  1983  under  the  direction  of  Dr.  R.  W.  Whalin,  Chief,  and 
Dr.  L.  E.  Link,  Jr.,  Assistant  Chief.  Dr.  L.  Z.  Hales,  Research  Hydraulic  En¬ 
gineer,  Mr.  K.  A.  Turner,  Computer  Specialist,  Ms.  M.  L.  Hampton,  Computer 
Technician,  Mr.  R.  E.  Ankeny,  Computer  Technician,  and  Mr.  K.  M.  Strausbaugh, 
Civil  Engineering  Technician,  performed  the  experimental  portion  of  the  study 
described  herein.  The  numerical  model  for  determining  wave  heights  downcoast 
of  a  straight  breakwater  at  an  angle  to  the  shoreline  under  combined  refrac¬ 
tion  and  diffraction,  based  on  the  uniformly  valid  asymptotic  theory,  was  ob¬ 
tained  by  contract  with  Dr.  Philip  Liu,  Cornell  University.  Dr.  Hales  pre¬ 
pared  this  report. 

Commanders  and  Directors  of  WES  during  the  conduct  of  this  investigation 
and  the  preparation  and  publication  of  this  report  were  COL  Nelson  P.  Conover, 
UK,  and  COL  Ti 1  ford  U .  (reel,  CE .  Technical  Director  was  Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  US  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


US  customary  units  of  measurement  used  in  this  report  can  be  converted  to 
metric  (SI)  units  as  follows: 
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feet 
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feet  per  second  per  second 


pounds-second-second  per  foot 
per  foot  per  foot  per  foot 
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0.3048 
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0.3048 


52.5540137 


0.09290304 
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second 
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EROSION  CONTROL  OF  SCOUR  DURING  CONSTRUCTION 


EXPERIMENTAL  MEASUREMENTS  OF  REFRACTION  AND  DIFFRACTION 
DOWNCOAST  OF  AN  OBLIQUE  BREAKWATER 


PART  I :  INTRODUCTION 


Statement  of  the  Problem 


1.  When  major  structures  are  erected  in  the  coastal  zone ,  they  alter 
currents  that  are  in  dynamic  equilibrium  with  the  existing  bathymetry.  These 
altered  currents  may  change  the  existing  bathymetry.  In  addition,  waves  break¬ 
ing  on  the  new  structure  will  cause  bottom  material  to  be  suspended  and  trans¬ 
ported  from  the  region  by  longshore  or  other  currents.  This  removal  of  mate¬ 
rial  from  around  structures  is  often  not  compensated  by  an  influx  of  additional 
material;  the  result  is  scour,  or  erosion,  that  usually  develops  along  the  toe 
of  the  structure.  In  order  to  ensure  structural  stability,  the  scour  area  must 
be  filled  with  nonerodible  material  (sufficiently  stable  to  withstand  the  en¬ 
vironmental  forces  to  which  it  will  be  subjected).  This  may  result  in  addi¬ 
tional  quantities  of  material  being  required  during  construction  that  can  po¬ 
tentially  be  very  costly.  To  minimize  potential  cost  increase  due  to  scour 
during  construction,  it  is  necessary  to  quantify  the  probability  and  ultimate 
extent  of  potential  scour  during  the  scheduled  construction  period. 

2.  Effective,  comprehensive,  and  low  cost  procedures  do  not  exist  for 
eliminating  scour  during  construction  in  the  nearshore  environment.  Determina¬ 
tion  of  potential  alternative  procedures  is  seriously  hampered  by  the  inability 
to  predict  the  extent  of  potential  scour.  Objectives  of  the  Erosion  Control 

of  Scour  During  Construction  research  program  are  to  develop  techniques  to 
minimize  and  control  scour  during  nearshore  construction,  and  to  predict  the 
probable  magnitude  of  scour  that  may  result  as  a  function  of  currents  and  wave 
climate.  One  phase  of  the  research  program  is  development  of  numerical  tech¬ 
niques  (incorporating  both  refraction  and  diffraction  effects)  for  computing 
the  wave  field  in  the  vicinity  of  structures. 

3.  The  present  state  of  nearshore  current  and  wave  theories  has  reached 
the  point  where  detailed  experimental  investigations  are  required  for  the 
verification  of  analytical  developments  and  numerical  models.  To  provide  a 


4 


firm  foundation  for  further  advancements,  a  simple  beach  profile  consisting  of 
straight,  uniform  contours  parallel  to  the  shoreline  was  experimentally 
studied  by  Hales  (1980).  A  shore-connected ,  vertical,  thin,  impermeable  bar¬ 
rier  (breakwater)  was  installed  perpendicular  to  the  shoreline  to  simulate 
prototype  jetties  and  breakwaters  commonly  occurring  along  many  coasts.  This 
single  jetty  (shore-connected  breakwater)  simplified  the  experiment,  facili¬ 
tated  direct  comparisons  with  numerical  model  results,  and  provided  greater 
understanding  and  insight  into  the  phenomenon  of  wave-height  variations  down- 
wave  of  a  breakwater  than  would  a  more  complex  geometry. 

4.  Analytical  models  of  wave  fields  surrounding  shore-normal  break¬ 
waters  have  been  developed  based  on  asymptotic  theory  and  the  parabolic  ap¬ 
proximation  (for  example,  Liu  and  Mei  1975,  1976;  Liu,  Lozano,  and  Pantazaras 
1979;  and  Lozano  and  Liu  1980).  These  analytical  models  have  been  compared 
with  the  experimental  data  of  Hales  (1980)  by  Liu  (1982)  and  Tsay  and  Liu 
(1982),  and  the  overall  agreement  between  theory  and  experiment  was  considered 
to  be  good.  Knowledge  of  these  important  phenomena  can  be  used  as  the  basis 
for  advanced  studies  of  sediment  movement  around  major  structures  under  com¬ 
bined  effects  of  refraction  and  diffraction,  when  the  structure  is  oriented 
perpendicular  to  the  shoreline. 

Objective  and  S cop e  of  the  Study 

5.  The  objective  of  the  present  study  was  to  extend  the  previous  experi¬ 
mental  work  of  Hales  (1980)  by  installing  a  shore-connected,  vertical,  thin, 
impermeable  barrier  (breakwater)  at  a  60-deg  angle  to  the  shoreline  to  simu¬ 
late  a  larger  range  of  prototype  jetties  and  breakwaters  in  existence  and  to 
test  a  new  theory  of  waves  in  the  lee  of  a  structure.  The  experimental  mea¬ 
surements  were  made  to  obtain  quantitative  knowledge  of  combined  refraction 
and  diffraction  in  the  lee  of  a  jetty  or  shore-connected  breakwater.  These 
data  were  then  compared  with  the  uniformly  valid  asymptotic  theory  of  Liu, 
Lozano,  and  Pantazaras  (1979)  for  the  same  arrangement,  thus  providing  verifi¬ 
cation  data  for  this  numerical  approximation  which  has  not  been  previously 
available.  This  report  presents  details  of  the  experimental  investigation  to 
measure  combined  refraction  and  diffraction  in  the  lee  of  a  jetty  or  break¬ 
water.  A  literature  review  of  the  theory  of  refraction  and  diffraction  is 
presented.  In  addition,  a  new  uniformly  valid  asymptotic  theory 


of  combined 


PART  II:  EXPERIMENTAL  INVESTIGATION 


Exper imen tal  Layout 

6.  This  study  to  investigate  wave  heights  downcoast  of  an  obliquely 
oriented  shore-connected  breakwater  or  jetty  under  the  combined  effects  of 
refraction  and  diffraction  was  conducted  in  the  experimental  facilities  of 
the  US  Army  Engineer  Waterways  Experiment  Station  (WES).  The  experimental 
arrangement,  which  was  molded  in  cement  mortar,  consisted  of  a  50-  by  60-ft 
area,  with  a  water  depth  of  1  ft  in  the  open-ocean  region  (Figure  1).  The 


"  A  table  of  factors  for  converting  US  customary  units  of  measurements  to 
metric  (SI)  units  is  presented  on  page  3. 


9.  The  ADACS  are  capable  o 1  automatically  calibrating  the  wave  sensors, 
controlling  wave  generators,  acquiring  data  from  the  sensors  at  a  high  sampling 
rate,  and  analyzing  test  data.  Data  are  taken  and  recorded  on  disc  or  magnetic 
tape  for  direct  analyses  by  the  minicomputer  system  or  on  magnetic  tape  in  a 
format  compatible  with  a  Honeywell  DPS1  for  backup  analyses.  Automatic  cali¬ 
bration  of  wave  sensors  has  reduced  the  time  required  to  calibrate  the  sensors 
by  a  factor  of  four.  In  addition,  several  times  the  number  of  tests  can  be 
run  during  a  day  with  test  results  analyzed  at  completion  of  model  tests  by 
minicomputer.  The  system  configuration  (Figure  2)  of  ADACS  consists  of  the 
following  subsystems:  (a)  digital  data  recording  and  controls;  (b)  analog  re¬ 
corders  and  channel  selection  circuits;  (c)  wave  sensors  and  interfacing  equip¬ 
ment;  and  (d)  wave  generators  and  control  equipment. 

Wave  sensors 

10.  The  data  acquired  from  wave  experiments  are  the  water-surface 
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Figure  2.  Schematic  of  components  of  ADACS 


each  set  of  parallel  roils,  the  voltage  from  the  signal  conditioning  equipment 
is  monitored  and  recorded  as  the  parallel  rods  are  moved  vertically  a  known 
distance  into  or  out  of  the  water.  A  precision,  linear-position  potentiometer 
is  located  on  the  wave  sensor  stand  and  is  coupled  directly  to  the  parallel 
rods  by  a  gear  train  driven  by  an  electric  motor.  By  moving  vertically  the 
coupled  wave  sensor „ind  potentiometer  with  the  electric  motor  and  by  monitor¬ 
ing  the  output  voltage  from  the  potentiometer,  the  wave  sensor  can  be  moved 
vertically  a  precise  distance.  The  electric  motor  for  each  wave  sensor  is  con¬ 
trolled  by  a  control/sense  line  and  a  relay  contact.  The  minicomputer  controls 
the  vertical  movement  of  each  wave  sensor  by  actuating  the  control /sense  line. 
The  central  processing  unit  acts  as  a  voltage  comparator  by  monitoring  the  po¬ 
tentiometer  voltage  and  comparing  it  with  a  reference  voltage  which  is  deter¬ 
mined  from  desired  displacement  and  potentiometer  calibration.  When  the 
voltage  comparison  is  satisfied,  the  control/sense  line  is  reactivated,  the 
electric  motor  stops,  and  voltage  samples  from  the  rods  and  potentiometers  are 
acquired.  By  systematically  moving  each  wave  sensor  through  11  quas i -equa 1 1 y 
spaced  locations  over  the  range  of  rod  length  used,  voltage  versus  known  dis¬ 
placements  are  obtained  from  which  a  calibration  curve  for  each  sensor  can  be 
calculated  and  recorded  on  magnetic  tape  or  disc.  After  collecting  the  cali¬ 
bration  data,  the  minicomputer  analyzes  these  data  by  least-squares  fitting  a 
set  of  curves  (linear,  quadratic,  or  spline)  to  the  data,  determining  the  best 
order  of  fit,  and  comparing  the  maximum  deviation  of  the  best  fit  with  a  pre¬ 
viously  acceptable  value  for  this  maximum  deviation. 

Data  acqu i s i t i on  and  a n a  1  vs i s 

12.  During  the  acquisition  mode,  wave  data  for  a  specified  wave  condi¬ 
tion  at  the  wave  generator  are  collected  from  a  maximum  of  50  wave  sensors, 
recorded  on  analog  strip  charts,  digitized,  and  recorded  on  magnetic  tape  or 
disc  for  further  analyses.  The  sampling  scheme  is  flexible  and  can  be  tailored 
for  different  applications  with  maximum  throughput  rates  theoretically  limited 
by  the  multiplexer  rate  and  allocatable  buffer  size.  The  sampling  scheme  used 
in  this  investigation  was  60  discrete  voltage  samples  equally  spaced  over  each 
wave  period  for  a  predetermined  number  of  90  wave  periods  for  each  of  the  sen¬ 
sor  locations.  The  minicomputer  calculated  from  input  parameters  the  lag  at 
the  beginning  of  data  acquisition  by  10  wave  periods  after  starting  the  gen¬ 
erator,  provided  timing  pulses  for  synchronizing  and  controlling  the  recorders, 
and  determined  completion  of  the  test.  The  determination  of  the  height  of 
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each  wave  of  the  monochromatic  wave  train  was  performed  (at  each  sensor  loca¬ 
tion),  the  average  of  these  90  individual  heights  was  calculated,  and  the 
standard  deviation  of  these  individual  observations  about  the  mean  was  com¬ 
puted.  The  value  ultimately  determined  as  the  wave  height  at  each  sensor  lo¬ 
cation  was  this  mean  value  plus  or  minus  one  standard  deviation 

Wave  Gage  Locations 

13.  The  wave-height  data  downcoast  of  the  experimental  breakwater  th.it 

was  positioned  at  a  6C-deg  angle  with  the  shoreline  were  obtained  by  operating 
36  wave-height  sensors  simultaneously  for  each  test  condition.  Two  ol  these 
wave-height  sensors  were  located  in  the  deeper  water  I  1  It  deep)  near  the  wave 
generator  to  ascertain  the  initial  generated  wave  height.  lhe  remaining  f4 
sensors  were  positioned  along  four  lines  parallel  with  the  shoreline  at  dis¬ 
tances  of  6,  8,  10,  and  13  ft  f rom  the  shoreline  (Figure  4);  the  still-water 

depths  at  these  four  sections  were*  0.3,  0.4.  0.3,  and  0.6  ft,  respectively. 

The  wave  gages  were  placed  on  a  supporting  pi  it  form  such  that  only  the  wave 
sensor  probe  penetrated  the  water  surface ,  thereby  eliminating  any  local  dis¬ 
turbance  caused  by  instrument  stands  touching  tin1  water  surface.  A  dry  bed 
photograph  of  the  experimental  facility  with  the  experimental  breakwater  in 
place  is  shown  in  Figure  3. 

14.  In  order  to  compare  experimental  results  with  numerical  or  analyti¬ 
cal  investigations,  it  is  necessary  t<>  have  goad  definition  of  the  experimental 
data.  The  34  wave-height  sensors  were  initially  placed  at  2-it  intervals  along 
the  four  sections  parallel  w i t h  the  shoreline,  and  the  complete  set  of  experi¬ 
mental  data  was  obtained  (gage  arrangement  No.  I,  Figure  0  )  .  To  provide  a 
more  dense  data  display,  the  entire  physical  arrangement  of  the  wave  sensors 
was  displaced  laterally  along  the  four  section  lines  tor  a  distance  of  1  ft, 
and  the  same  set  of  wave  data  (period  and  height)  was  repeated.  Ihis  second 
testing  arrangement  is  shown  as  gage  arrangement  No.  3,  Figure  7.  All  data 
from  these  two  gage1  arrangements  resulted  in  a  data  set  that  presented  the 
wave  heights  at  1-ft  increments  along  the  four  sections  which  run  parallel 
with  the  shoreline,  with  the  data  extending  from  very  near  the  breakwater  deep 
in  the  shadow  zone  and  extending  across  the  region  where  the  waves  experience 
both  refraction  and  diffraction.  These  precise  experiments  provide  data  that 
define  the  wave  height  downcoast  of  the  structure  and  can  he  used  for  compari¬ 
son  with  numerical  or  analytical  studies  of  the  same  concept. 

13 
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PART  III:  LITERATURE  REVIEW  OF  COMBINED 
REFRACTION  AND  DIFFRACTION 


15.  The  approximate  solution  of  water  wave  refraction  caused  by  a  vari¬ 
able  bathymetry  is  well  known  and  can  be  derived  by  assuming  that  bottom  re¬ 
flections  are  negligible  or  by  the  more  rigorous  Wentzel -Kramers-Br i 1 louin 
(WKB)  approximation.  The  exact  solution  for  the  diffraction  of  surface  waves 
by  vertical  barriers  of  simple  cross  section  in  water  of  constant  depth  also 
is  well  known  (being  analogous  to  classical  problems  of  physics).  However,  an 
analytical  theory  for  the  practical  case  of  combined  refraction  and  diffrac¬ 
tion  has  not  been  completely  developed.  Th-  present  engineering  practice  for 
determining  wave  heights  under  this  condition  is  a  stepwise  procedure  (Dunham 
1951,  Liu  and  Mei  1976,  US  Army  Coastal  Engineering  Research  Center  1977). 

The  procedure  involves  the  following  steps:  (a)  calculate  refraction  effects 
up  to  the  barrier  (jetty  or  shore-connected  breakwater);  (b)  calculate  for  a 
"few"  wavelengths  the  effects  of  diffraction,  assuming  a  constant  water  depth; 
and  (c)  beyond  this  region  calculate  the  refraction  effects  only.  This  pro¬ 
cedure  is  obviously  imprecise,  and  Mobarek  (1962)  indicates  that  the  method  is 
suitable  only  for  intermediate  water  depths.  In  addition,  Whalin  (1972) 
states  that  this  procedure  is  only  valid  for  small  refraction  effects. 

Wa t e r_Wa ve  Ref rac ti on 

16.  In  intermediate  and  shallow  water,  the  phase  speed  of  a  surface 
gravity  wave  depends  on  water  depth.  Since  wave  celerity  decreases  as  depth 
decreases,  phase  velocity  varies  along  the  crest  of  a  wave  propagating  at  an 
angle  to  underwater  contours  because  that  part  of  the  wave  in  deeper  water 
moves  faster  than  that  part  in  shallower  water.  This  variation  causes  the 
wave  crest  to  bend  toward  alignment  with  the  contours.  This  bending  effect 
due  to  changes  in  bottom  topography,  called  refraction,  depends  on  the  rela¬ 
tion  of  water  depth  to  wavelength,  d/L  ,  and  is  analogous  to  refraction  of 
other  types  of  waves  such  as  light.  A  basic  assumption  in  wave  refraction 
theory  is  the  conservation  of  energy  between  wave  orthogonals  (i.e.,  no  dif¬ 
fraction  of  energy  along  wave  crests).  The  change  in  wave  direction  of  dif¬ 
ferent  parts  of  the  wave  results  in  convergence  or  divergence  of  wave  energy 
and  materially  affects  the  forces  exerted  by  waves  on  structures  and  of  the 


capacity  of  waves  t  «>  l  i  ansport  sand  eiltui  alongshore  or  oiisln.i  r/nl  I  shore. 

17.  I’rcu  (■•in  res  ter  !  lit*  compu  t  a  t  i  mi  el  re ! r  a .  t  i on  •  >  t  mii  lair  gravity 
waves  un  water  el  neminiterm  depth  involve  tin  i : .  sump  t  i  on  that  a  wave  with  a 
curved  f  rest  pattern  ami  variable  amp  1  i  t  lute  li.uig  !  lie  i  list  behaves  locally  as 
a  st ra i ght -c rest ed  wave  oi  constant  amplitu.le.  Kivleigh  (  1 8  7  7  )  appears  to 

have  Inn  ui  the  I  i  rst  to  use  the  a  pp  rox  I  m  1 1  i  •  ui  s  I  . . .  t  n'al  apt  us  in  this 

analysis,  amt  theoretical  results  have  t-eeu  lev,  l-ped  with  lospeit  to  energy 
flux  ami  phase  speed.  As  expressed  hv  Keller  I  i 11  i  .  tin  geometrical  opt  its 
theory  defines  a  propagation  ve  lot  I  t  y  at  each  point  ui  tin  water  ini  late,  with 
this  velocity  being  exact  Iv  thit  which  w.tv,*s  o)  giv-n  pi-iiod  would  have  ill 
water  id  uniform  depth  at  ail  points.  Hv  employing  Fermat's  principle  ol 
optics,  wave  ra Vs  a i e  del i net  and  suria-e  waves  are  assumed  to  propagate  along 
these  rays.  The  variation  of  t  tie  amplitude  a  long  the  rays  is  determined  hy 
the  use  of  tin  prim  tple  ohm  aval  mn  of  energy.  Ibis  principle  (in  its 

•  •pt  i  ra  1  formi  slate;  that  l  f  1  tix  •>!  energy  is  the  same  at  all  cross  sections 
between  two  adjacent  warn-  rats.  i  In*  energy  I  lax  is  proportional  to  the  square 
ot  the  amplitude  d  tin  wives  and  to  t  lie  distant,  between  the  rays,  and  hence 
the  wave-he  i  gilt  va  r  i  a  t  i  n  along  the  i  a  v  is  a*.  a  i  1  aide  . 

18.  Keller’s  (  I  h  TK  i  lei  iv.it  mn  is  ha  so  1  upon  an  "asymptotic"  solution 
o  |  the  equut  ions  of  the  ex  a  i  t  linear  l  h>  :  y  f  ( •  i  p,  i  i  ■  "i  i  t  waves  in  wafer  •>! 

a  r  h  i  t  ra  rv  nonun  i  1  ■■  rin  •  n  p  t  n  i  lie  vdul  i  n  ■  s  .  vinp  t  •  ■  t  i  •  in  flu-  seme  that  the 


depth  and  wavelength  a  i  <  .in  i 


t  red  with  t  he  hoi  i  .•  nt  a  1  s,  si  1 1  f  the  hot  torn 


contours  .  l'he  first  term  ■ ,  I  i  I  i.l  i  ,.n  atm"  t  i  v  with  tin  i  s  \  nipt  o  t  i  c 

t  o  rm  of  t  fie  solution  tor  waves  in  w  a  t  "  r  „  i  t !;  a  un  i  I  rm  1  y  sloping  ho  t  t  mn  as  the 

bottom  slope  tends  to  zero.  i  h  i  s  solution  ,  on  I  i  ms  with  .ill  the  principles  of 
t  he  .geonie  t  r  i  c  ,i  1  optics  theory  of  Kay  1  <•  i  gh  <  IK/'/)  ml  tints  provides  a  deriva¬ 
tion  of  that  t  heo rv  lie-  results  are  not  valid,  however,  at  c  a  u  x t  i  ,  s  or  ray 

in  oss i ngs  . 

1  ‘1 .  In  pro!.  |,  iic  ,,t  .in*. ii  wave  propag.it  •  •  ■  1 1  over  mild  slopes,  the  prin- 

i  i  p  I  e  of  geomo  t  r  :  ■  a  1  his  been  applied  by  (artier  (  1  Ohb )  as  the  1  i  rs  t 


a  pp  rox 


liiiation  in  a  systematic  perturbation  scheme  while  the  bottom  is  con¬ 


sidered  to  fie  i  oe  a  I  1  v  tu ,  i  i  /o  u  t  a  I  .  1  he  depth  variation  was  dealt  with  after¬ 
ward  hy  requiring  the  appropriate  energy  conservation.  Ibis  was  a  1  so  the 
basic  idea  for  tin*  work  of  Kofi  and  heMc  haute  (  I  dfif,  i  in  which  the  transform.'!- 
lion  of  progressive  waves  wts  investigated  as  t  lie  c  travel  from  deep  water  to 
shore.  The  Stokes'  l  heo  i  v  at  a  fifth  "pier  <  ■  ‘  a  pp  i  •  ■  <  i  n,a  t  t  "'!  was  applied  a  I  on; 


with  the  method  of  cmiservat ion  of  energy  flux.  It  is  assumed  that  the  wave 
is  simply  harmonic  in  time.  The  first,  third,  and  fifth  orders  of  approxima¬ 
tion  were  compared  with  each  other  and  with  experiments.  The  differences  be¬ 
tween  the  predictions  of  wave-height  changes  based  on  the  three  orders  of  ap¬ 
proximation  were  found  to  he  small,  on  the  order  of  5  percent.  For  practical 
purposes,  the  third-order  theory  was  found  to  give  reliable  results.  The 

third  and  fifth  order  Stokes'  theories  are  based  on  a  series  expansion  in 

3  5 

terms  of  H/  I.  where  terms  of  t  fie  order  of  (H/L)  and  (H/L)  ,  respectively 

are  retained  and  higher  order  terms  are  neglected.  it  should  be  noted  that 
this  theory  is  based  on  an  expansion  in  term  of  the  wave  steepness,  H/L  ,  and 
consequently  can  he  expected  to  better  approximate  limit.ag  steepness  waves  in 
deep  water.  However,  it  cannot  be  expected  to  adequately  predict  wave  cuarac- 
teristics  in  shallow  water,  since  water  depth  is  not  a  parameter  in  the  series 
expansion.  Thus  the  theory  is  a  finite  amplitude  deepwater  wave  theory. 

20.  In  cases  of  limiting  shallow  water,  the  wave  conditions  are  nearlv 

independent  of  wavelength,  and  the  important  parameters  are  water  depth  and 

the  ratio  of  wave  height  to  water  depth.  Keulegan  (1950)  showed  that  Stokes 

waves  are  most  nearly  valid  in  water  deeper  than  about  d/L  >  1/8  to  1/10  . 

In  shallower  water,  cnoidal  wave  theory  appears  to  be  more  satisfactory,  and 

Masch  (1964)  investigated  the  problem  of  wave  shoaling  using  cnoidal  wave 

theory  with  the  formulas  developed  by  Keulegan  and  Patterson  (1940).  Masch 

(1964)  assumed  hydrostatic  pressure  distribution  and  neglected  the  convective 

inertia  term  in  his  expression  for  the  energy  flux.  The  third  and  fifth  order 

cnoidal  theories  are  based  on  a  series  expansion  in  terms  of  H/d  ,  where 

3  5 

terms  of  the  order  of  (H/d)  and  (H/d)  ,  respectively,  are  retained  and 

higher  order  terms  neglected.  It  should  be  noted  that  this  theory  is  based  on 
an  expansion  of  the  relative  wave  height  (H/d)  and  can  be  expected  to  better 
approximate  the  wave'  form  in  shallow  water.  However,  it  cannot  fie  expected 
to  do  a  very  good  job  ot  approximating  the  wave  form  for  limiting  steepness 
waves  in  deep  water.  in  that  case,  water  depth  is  unimportant  and  wavelength 
is  crucial.  This  tfieory  could  he  considered  a  finite  amplitude  intermediate 
and  shallow-water  theo.y. 

21.  A  technique  of  "asymptotic  expansion"  was  applied  by  Mei,  Tlapa, 
and  Eagleson  (1968)  to  water  wave  propagation  over  an  uneven  bottom  that  has 
straight  and  parallel  confiiirs.  Attention  was  focused  on  the  establishment  of 
a  rigorous  scheme  of  successive  approximation  tor  higher  order  corrections. 


g  =  gravitational  constant,  32.174  ft/sec 
k  -  wave  number,  2n/L  ,  1/ft 
d  =  local  water  depth,  ft 
a  =  local  wave  amplitude,  ft 

The  second  derivative  of  the  wave  amplitude  in  the  horizontal  plane  is  given  bv 
a^^  or  a  .  The  rate  of  power  transmission  P  ,  or  energy  flux,  was  deter¬ 
mined  to  be: 


=  G pga2n  k Ab)  f1  + 


where 

Ab  =  wave  ray  spacing,  ft 

2  4 

p  =  fluid  density,  lb-sec  /ft 

n  =  ratio  of  group  velocity  c^  ,  to  phase  ,'elocity  c  ,  dimensionless 
to  =  angular  frequency,  2rt/T  ,  1/sec 


Define : 


a  .  (^) 

The  commonly  used  existing  procedures  (Dunham  1951,  Liu  and  Mei  1976,  CERC 
1977)  for  construction  or  computation  of  refraction  diagrams  utilize  phase 
speeds  that  are  obtained  by  neglecting  6  . 

23.  Battjes  (1968)  examined  the  omission  of  6  from  a  qualitative 
standpoint.  In  an  area  of  strong  local  convergence,  omission  of  6  from  Equa¬ 
tion  1  results  in  underestimating  the  local  phase  speed.  The  result  is  that  in 
Equation  2  the  ray  separation  Ab  will  be  underestimated.  Thus  there  are  gen¬ 
erally  two  contributions  to  the  error  that  results  in  the  computed  wave  ampli¬ 
tude  a  .  However,  these  two  contributions  are  of  a  different  nature  because 
the  effect  on  the  wave  pattern  of  using  an  incorrect  value  for  the  phase  speed 
is  cumulative,  whereas  the  effect  on  energy  flux  is  local.  In  any  case,  the 
omission  of  6  will  generate  wave  heights  at  variance  with  the  height  inferred 
from  refraction  diagrams  based  on  linear  small  amplitude  wave  theory  which  ne¬ 
glects  the  effects  of  wave-height  gradients  along  the  wave  crest,  so  that  it 
would  appear  that  energy  had  been  transferred  across  orthogonals.  An  estimate 
of  the  magnitude  of  6  has  been  approximated  by  Battjes  (1968)  for  four 


distinctly  different  cases.  For  simple  shoaling,  the  amplitude  varies  only  in 
the  direction  at  wave  propagation,  and  tor  shallow  water  the  variation  was: 


6  -  8  * 


14) 


where 

L  =  local  wavelengths 

s  =  bottom  slope 

For  short-crested  waves  where  the  distance  along  the  wave  crest  is  two  or  three 
times  the  wavelength,  6  =  -20  percent  or  -10  percent,  respectively.  For  the 
case  of  dif  fract  ion  around  a  serai -ini  ini te  breakwater,  6  was  found  to  reach 
values  between  +10  percent  and  -10  percent  at  points  a  distance  of  one  wave¬ 
length  from  the  breakwater  tip,  decreasing  inversely  proportional  to  the  dis¬ 
tance  from  the  tip.  For  the  case  of  diffraction  through  a  gap  of  width  two 
times  the  wavelength,  6  was  evaluated  in  a  few  points  on  the  center  line  of 
the  gap  where  it  was  found  to  reach  values  of  25,  7,  and  3  percent  at  distances 
of  1,  2,  and  4  wavelengths,  respectively,  from  the  gap. 

Water  Wave  Pi f  f  rac t i on 

24.  Diffraction  of  water  waves  is  the  phenomenon  by  which  wave  energy 
propagates  into  the  sheltered  lee  of  structures  even  in  the  absence  of  bathy¬ 
metric  refraction.  in  these  situations,  wave  crests  bend  (even  in  constant 
depth  water)  and  gradients  of  wave  height  exist  along  the  wave  crest.  This 
phenomenon  is  most  visible  when  a  train  of  regular  waves  is  interrupted  by  an 
obstruction  such  as  a  jetty  or  shore-connected  breakwater.  The  theory  of  water 
wave  diffraction  can  be  explained  by  Huygens'  principle.  Each  point  of  an  ad¬ 
vancing  wave  front  (wave  crest)  may  be  considered  as  the  center  of  a  secondary 
circular  wave  which  advances  in  all  directions.  The  resultant  shape  of  the 
crest  is  the  envelope  of  all  these  secondary  waves.  In  a  straight-crested  wave 
train,  the  envelope  of  the  secondary  waves  is  a  straight  line  also.  When  the 
wave  passes  an  obstruction,  the  energy  intensity  at  a  certain  point  is  a  vector 
combination  of  all  the  circular  waves  emitted  by  every  point  of  the  passing 
wave  train. 

25.  Somme r f e 1 d  (1896)  presented  a  solution  for  the  diffraction  of  light 
waves  past  the  edge  of  a  semi - i n f ini t e  screen.  Penny  and  Price  (1944)  showed 
that  this  is  also  the  solution  of  the  water  wave  diffraction  problem  at  the 


Figure  8.  Definitive  sketch,  wave  diffraction 
around  a  breakwater 


end  of  a  semi - inf  ini te  obstacle  such  as  a  jetty  or  shore-connected  breakwater. 
This  exact  solution  of  the  surface  elevations  behind  the  breakwater  is  applica 
ble  only  to  water  of  constant  depth  and  waves  of  small  amplitude.  Putnam  and 
Arthur  (1948)  summarized  the  solution  of  Penny  and  Price  (1944)  for  the  de¬ 
finitive  sketch  of  Figure  8.  In  cylindrical  coordinates,  the  water-sur face 
elevation  is: 


/aikc\  ikct  ,  ,  , 

r|  =  ^  — — Je  cosh  kd  F(r,6) 


(5) 


where 

r|  =  water-surface  elevation,  ft 
i  =  square  root  of  -1 

and  the  other  symbols  have  been  previously  defined.  F(r,0)  is  a  function 
which  satisfies  the  wave  equation  in  cylindrical  coordinates: 


q  ,  1 1£ ,  i  §*  ,  a  ,  0 

:)r2  r  ar  r2  ae2 


(6) 


In  the  presence  ot  a  jetty  or  breakwater,  the  boundary  condition  is  imposed 
so  that  the  normal  component  of  the  fluid  velocity  is  zero  along  the 
breakwater,  leading  to  the  solution: 


1 

,  i((n/4)-kr  cos  (6  -6)]  f  2..,, 

F(r ,0)  =  e  °  /  e'(inu  /2)  du 

42  J 


,  i[(n/4)-kr  cos  (6  -0)) 

+  ±-  e 

42 


u , 

/ 


e'(i7Iu2/2)  do  (7) 


Bretschneider  (1966)  has  presented  computational  procedures  for  evaluating  the 
diffraction  coefficients  at  arbitrary  points  behind  jetties  or  breakwaters. 

26.  Wiegel  (1962)  developed  a  graphical  procedure  for  determining  dif¬ 
fraction  coefficients  of  waves  passing  the  tip  of  single  breakwaters.  The 
family  of  diagrams  shows,  for  uniform  water  depth,  lines  of  equal  wave-height 
reduction  displayed  in  terms  of  the  diffraction  coefficients.  The  diffraction 
diagrams  (typical  example,  Figure  9)  are  constructed  in  polar  coordinate  form 
centered  at  the  structure  tip.  The  arcs  behind  the  breakwater  are  spaced 

one  radius-wavelength  unit  apart  so  in  application,  a  specific  diagram  must  be 
scaled  up  or  down  so  that  the  particular  wavelength  corresponds  to  the  scale 
of  the  hydrographic  area  under  investigation.  The  set  of  diffraction  diagrams 
of  waves  passing  the  tip  of  a  single  breakwater  was  presented  by  CERC  (1977). 
Figure  9  is  the  configuration  of  the  wave  approach  direction  analogous  to  the 
physical  hydraulic  model  layout  used  in  this  experimental  investigation. 

27.  Laboratory  tests  were  performed  by  Harms  (1979)  to  investigate  the 
distribution  of  wave  heights  in  the  lee  of  a  breakwater  (shore-connected)  for 
waves  normally  incident  upon  the  structure  and  with  a  horizontal  bottom  both 
in  front  of  and  in  the  lee  of  the  structure.  In  general,  satisfactory  agree¬ 
ment  was  obtained  between  measurement  and  theory,  but  difftaction  theory  was 
not  found  to  be  conservative.  At  large  distances  in  the  shadow  zone,  measured 
wave  heights  consistently  exceeded  theoretical  values.  Close  to  the  break¬ 
water  outside  the  shadow  zone,  the  measured  maximum  wave  height  was  also  found 
to  be  larger  than  that  predicted  by  theory.  The  diffraction  behavior  appeared 


to  be  insensitive  to  the  intensity  of  wave  reflections  from  the  seaward  side 
of  the  breakwater. 

Combined  Water  Wave  Refraction  and  Diffraction 

28.  The  bathymetry  shoreward  of  a  jetty  or  breakwater  usually  is  not 

flat  or  even  uniform;  hence,  refraction  generally  occurs  in  addition  to  the 
diffraction  effects.  While  a  general  unified  analytical  approach  to  the 
simultaneous  solution  of  these  two  distinctly  different  phenomena  has  not  been 
entirely  developed,  considerable  insight  has  been  gained  through  the  theo¬ 
retical  work  of  Liu  and  Mei  (1975,  1976),  Lick  (1978),  Liu  and  Lozano  (1979), 
and  through  the  earlier  experimental  work  of  Mobarek  (1962).  The  procedure 
usually  followed  by  coastal  engineers  concerned  with  wave-height  variation  be¬ 
hind  jetties  or  breakwaters  is  to  construct  refraction  diagrams  shoreward  to 
the  structure,  then  construct  diffraction  diagrams  for  three  or  four  wave¬ 
lengths  shoreward  of  the  jetty,  and  finally  refract  the  last  wave  crest  on 
toward  the  shoreline.  This  procedure  is  schematized  in  Figure  10  where  the 
overall  refraction-diffraction  coefficient,  ^  ,  in  the  region  behind  the 

structure  is: 

Kr-d  =  KrKd  ^ 

where 

K  =  refraction  coefficient  at  the  structure,  dimensionless 
r 

Kj  =  diffraction  coefficient  on  last  wave  crest  behind  the  structure 
from  which  additional  refraction  computations  are  performed, 
d imens ionless 

bj  =  orthogonal  spacing  at  the  last  diffracted  wave  crest,  ft 
=  orthogonal  spacing  near  the  shore,  ft 

29.  Mobarek  (1962)  experimentally  investigated  the  effect  of  bottom 
slope  on  wave  diffraction  through  a  gap  in  a  breakwater  normal  to  the  incident 
wave  direction.  Also  investigated  was  the  effect  of  an  abrupt  increase  or  de¬ 
crease  in  the  water  depth  behind  the  breakwater.  The  theoretical  analysis  for 
the  comparison  of  experimental  results  followed  the  treatment  of  Penny  and 
Price  (1944)  restricted  to  the  case  of  normal  incidence  for  which  the  Sommer- 
feld  (1896)  solution  is  simplified  and  in  the  presence  of  a  horizontal  bottom. 
Two  fundamentally  different  basin  configurations  were  used  in  the  study.  The 


26 


Figure  10.  Graphical  procedure  for  determining  overall  refraction  and 
diffiaction  effects  past  a  breakwater  (after  CERC  1977) 

first  consisted  of  a  longitudinally  sloping  bottom  with  the  slope  commencing 
at  the  breakwater  and  extending  to  the  shoreline.  The  second  was  constructed 
with  a  flat  bottom  extending  beyond  the  breakwater  gap  but  sloping  laterally 
to  the  shoreline.  Taking  into  consideration  the  serious  limitations  of  the 
experimental  equipment  (very  small  model,  72  sq  ft),  the  investigation  led  to 
the  conclusion  that  the  procedure  usually  followed  for  estimating  wave  heights 
behind  jetties  or  shore-connected  breakwaters  (Figure  10)  was  sufficiently 
good  for  medium  period  waves;  but  in  the  case  of  long  waves,  the  effect  of  the 
shoaling  bottom  on  waves  should  be  taken  into  consideration.  Experiments  on  a 
larger  scale  were  highly  recommended. 

30.  For  a  long  shore-connected  breakwater  on  a  slowly  varying  bottom, 
an  asymptotic  theory  has  been  developed  by  Liu  and  Mei  (1976)  that  accounts 
for  the  combined  effects  of  refraction  and  Fresnel  diffraction  of  water  waves. 
However,  for  short  jetties  or  groins,  the  reflection  and  diffraction  effects 


h'.  1 11.1 1  i  on  !  1  I'nlmi's  to  t  hr  diffraction  Hr  I  miio  1  t  /  equa  t  l  mi  in  deep  m  i  ■  >n  t  u  '  - 
depth  water.  In  shallow  water,  t  lie  equation  reduces  to  the  lineal  i.ng-we.i 
eijuat  ion  . 

33.  Houston  (1980)  solved  Hquat  mu  11  by  the  use  ol  a  hybrid  finite  .  le¬ 
nient  numencaL  model  originally  developed  by  Chen  and  Me  i  (1974)  to  solve  the 
diffraction  Helmholtz  equation  in  a  cons t ant -depth  region.  The  appropriate 
modifications,  including  variable  depth  and  frequency  dispersion,  were  in¬ 
corporated  by  Houston  11980)  and  the  solution  of  F.quation  11  was  applied  t, 

the  geometry  of  the  experimental  study  of  Hales  (1980)  (i.e.,  a  uniform  s|.|e 
with  a  shore-connected  breakwater  perpendicular  to  the  shoreline).  rublem 

in  simulating  those  experimental  hydraulic  tests  numerically  was  that  the 
waves  broke  m  the  experimental  facility  near  the  shoreline  and  thus  dissip.t.d 
their  energy.  No  mechanism  existed  to  dissipate  energy  in  the  numerical  m- ..  1- ■  1  - 
However,  dissipation  was  simulated  by  allowing  waves  to  continue  to  prepay,  ite 
out  of  tht1  problem  area.  I  he  breakwater  and  uniform  slope  were  simulated  ■ 
to  the  point  where  breaking  occurred.  The  depth  was  then  increased  to  the 
depth  of  a  semi - i n t i n i t e  ocean  region  surrounding  the  region  ot  cumpul at i  ■  n , 
and  the  waves  were  a  1 1  owed  to  radiate  away  from  the  area  of  interest. 

34.  Figure  11  shows  a  typical  comparison  between  the  experimental 
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figure  II.  Com  pa  i  i  si  ii  o  I  .isvinpt  nt  i  i  and  finite 

element  solutions  with  e:<pei  inieiila]  measurements 
(breukuulet  normal  to  shoreline)  (after  Houston 
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I  N  I  FuKMi.V  VALID  ASYMPTOTIC  THEORY 


Background 

ir>  In  11‘icnl  years,  the  parabolic  approximation  has  been  developed  ex¬ 
tensively  tor  studying  wave  scattering  problems  in  different  branches  of  the 
physical  sciences.  Hidde.  ,1979)  and  Lozano  and  Liu  (1980)  derived  indepen¬ 
dently  the  parabolic  approximation  for  water-wave  problems.  Analytical  solu¬ 
tions  were  obtained  f.-i  the  combined  refraction  and  diffraction  wave  field 
near  a  thin  breakwater  perpendicular  to  the  shoreline  on  a  plane  beach.  The 
background  wave  field  was  assumed  t a  have  straight-line  wave  rays. 

lb.  Also  based  on  the  parabolic  approximation,  a  numerical  study  of 
water-wave  retraction  and  diffract  ion  problems  has  been  conducted  by  Tsny  and 
Liu  (1982),  where  the  refraction  index  is  not  constant.  Two  problems  were 
considered:  (a)  the  wave  field  near  a  submerged  shoal  on  a  sloping  bottom  and 

(b)  the  wave  field  in  the  neighborhood  of  a  breakwater  on  a  sloping  beach.  In 
the  latter  problem,  the  orientation  of  the  breakwater  is  no  longer  limited  to 
be  perpendicular  to  the  shoreline.  For  the  perpendicular  case,  the  accuracy 
of  the  parabolic  approximation  numerical  results  was  verified  by  comparing 
with  the  precise  experimental  data  of  Pantazaras  (1979)  and  Hales  (1980)  (both 
sets  of  data  having  been  obtained  at  WES).  For  this  case,  the  uniformly  valid 
asymptotic  theory  of  Liu,  Lozano,  and  Pantazaras  (1979)  was  also  used  to  verify 
the  parabolic  approximation. 


Deve 1 opment  of 


37.  For  the  case  of  a  shore-connected  breakwater  on  a  linear  plane 
beach,  Liu  and  Me l  (197b)  and  Lozano  and  Liu  (1980)  have  previously  shown  that 
an  approximate  closed  form  solution  can  be  obtained  by  the  parabolic  approxi¬ 
mation.  However,  Liu  (1982)  showed  that  this  solution  becomes  invalid  near 
the  tip  of  the  breakwater  because  of  the  inherent  nature  of  the  parabolic  ap¬ 
proximation.  To  remove  this  weakness,  Liu,  Lozano,  and  Pantazaras  (1979)  and 
Liu  (1982)  developed  a  a  uniformly  valid  asymptotic  solution  for  the  same 
problem.  This  theory  is  verified  by  the  experimental  data  of  the  present 
study.  The  beach  topography,  which  is  required  to  be  uniform  in  the  along¬ 
shore  direction,  can  he  arbitrary  in  t  fie  onshore-offshore  direction. 


' 


P 


I 


i  • 


17 1.’  »■>' 


.58 .  Following  the  development  of  Liu,  Lozano,  and.  Pant  iz  me  1  I')/"1,  ate 
m  the  notation  at  Liu  (1982),  the  geometry  of  the  hreakwa!  er  i :  genera :  «/**d 
to  he  one  of  radiated  wave  rays  emitted  from  the  t 1 p  id  the  *m  <  ikwatm  (fig~ 

11 1  e  12).  The  relation  between  the  incident  wave  angle  and  (ho  ret  !e.  t  wave 
angle  along  the  breakwater  is  shown  in  Figure  Id.  Sma 1 l-.mip 1  . t  ud- 

S  r*  r'r 


!  t '  M  t 


Br  eQKwote< 


\ 


\  B  r  *  2  k  m  0  ’  f  1 


Reflected  wove 


Shoreline 

V 


1-22 

"  l  -  i — 
h  {  I  ) 

"'a 


R q d  -  j  t  e d  «o»f  r  2 ' 


incident 
wcwe  rOf 


\ 

e:\ 


figure  12.  Definitive  sketch, 
ini .  f  da  r'm  1  y  va  lid  1  symp  tot  ic 
t  he,  rv  1 1  i n ,  Lozano ,  and 
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Figure  13.  Definitive  rket'h.  uni 
formly  valid  asym;totic  t  h".a  v, 
relation  between  lie  idem  wave- 
angle  and  the  n  fin  ten  wave  angle 
along  the  breakwater  i  I ..  .  i.*»/.tn«#  , 

Fanla/a 1  as  1 9  1  • 
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wives  with  the  Me  I  dent  wave  amplitude,  a  ,  and  radian  ti’nien  v 

o 

,11. Slimed  in  the  development.  The  angle  ol  incident*'  is  del  in-'  1  is 
l.u/ario ,  and  Pant  izaras  ( .1979)  have  shown  that  the  1-ading  "f  l<  1  •>vt:ipt  el  1 
so  1  ut  ion  lor  the  t  re*  surface  displacement,  f)  ,  is. 
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A(x)  represents  the  combined  refraction  and  shoaling  factor: 
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The  subscript  o  in  Equation  lb  denotes  the  quantities  associated  with 
dent  waves  in  the  far  field.  The  function  G ( p ^ )  defined  in  Equation 
be  given  in  terms  of  the  Fresnel  integrals,  whence: 


(16) 

inci  - 
15  can 

(17) 


(18) 


(19) 


are  the  Fresnel  cosine  and  sine  integrals,  respectively.  The  arguments  of  the 
function  G(p  J  in  Equation  15  were  defined  by  Liu  (1982)  as: 
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function  (S  ,  S  ,  >r  K)  it  .mv  ,nin  t  r.n  v  j  ■>  <  >  i  n  t  can  be  tuns  i  dered  .is  the  sum 
of  the  wave  number  comp. uniit.  in  tfie  radial  direction  between  the  irbitrary 
point  and  Lin'  t  i  p  «.f  the  bre.ikw.it  er .  'I  lie  va !  in  ut  the  ph  ise  timet  ion  for  radi¬ 
ated  waves  at  any  urletrarv  point  it.  always  gt  eat  ot  than  or  eipia  i  to  S  and 
S  .  The  brnnehes  of  t  fie  mu  1  t  j  va  I  tied  (muttons  wi  and  U  are  del  tried  by  Liu 
(1982)  as  follow:-:  I  he  value  ot  U  is  ivy  ativ.  .ns.de  t  lie  shadow  region  de¬ 

fined  according  to  the  (*pnmct  r  i  c.i  1  opt  <  <.  s  I  ■  ry  and  is  positive  elsewhere. 


.here 


d  M  [ 


4 ( )  .  I  tlf 


1  u 

t  he  r 

r  f  :  tM.  1  j 

t  ri  rr g  l  uti 

1  fid 

is  n<-ga  t  i  .'<•  <  *  i  sew  In 

■  re  .  In  t.  lie 

h  \~t 

■a  kw  a  t 

r  r  hr. 

.1  fill  Vft! 

'i-'il  ('•' 

a!:d  I'Oilll'lde..  wifi 

-  one  of  tile 

r.i  y 1 

!  rom 

t  h.r  1  l 

!  [ 1  ,  l hr  ph  i 

i  *  1 

:  i.ro  t  i  on  for  t  he  re  1 

i.'i  t  ed  wave  i  s  : 

r 

X 

■! 

,  ( '  . 

!  .\  c  ii  v 

i  ' 

so) 

K  -  k 

i  i 

*.  i  »i 

{  1 

1  22  ) 

*  ■  1 

•:lt 

ilr  ur 

f  r..i!iif 'r!'  . 

•  t  -■< !  .i  ■  ■  •  o . i:  i  ■ 

i  !  'li,'  t  '|.' 

SOrM 

I'li 

1  1  ■  e  1 

v  . 

« i  i  •  r  •  \  t 

1  I  t.p  ij# ■  n  t  i  f 

(  i  • ' 

•  *  i  •  ••  :  i  "  1  .  ■  . 

. ; .  .  .!;•  1 1 1  doWM- 

n  “« 

o  n  n  >  ■  i 

If!  t:|. 

-.1  kv  :  t.  f  ■  .  ! 

1 .  j  :  >  f  • 

1  ■  ri  (  1. •  on  i  !  ■  :  el :  \ 

•  i  1  ; d  i  s vmpt ot  i ( 

<n 

12), 

.is  h.i:o 

•d  on  t  hr  ! 

1  l  f !  i  i 

Illie’i  t  i  ]  .i  ,  MHiip  l  '  !■  t 

hit  the  !  i  l  e  i  i:  - 

1  i : 

plain 

■  r  ro  •  i 

i  !  i  i  W  r  <  1  .  -. 

,  1  i  .it 

.  .  1  Wive  .  ]  \  1 .  1  '  '  > ' ' 

i  :  ;  ■  .»>  I  be  tip 

:  r  } 

J  J ,  f 

•  if  !  ,i  ; 

t  d  tiff  1  v  l  1 

it-li 

•  1  1  ti  1  :  .1 1 ;  I  i  ' '  1  ' 

■  ■  r '  i  "us 

required  in  the  deve  1  opnienl  .ire  presented  in  Liu,  Lozano,  and  Pantazaras  (  l(J7d) 
However,  m  actuality,  tf  ■  plan  form  layout  of  most  sho re- Conner t ed  breakwaters 
at  an  angle  to  the  shoreline  is  that  of  a  straight  line.  Hence,  while  the 
assumption  ot  a  planform  layout  following  a  wave  ray  radiated  tiom  the  tip  ot  tin 
structure  expedited  t  tie  uiiulytiiul  development,  the  numerical  model  based  on 
this  development  does  not  precisely  conform  with  most  prototype  conditions.  It 
was  therefore  desirable  to  adapt  the  numerical  scheme  to  fit  the  case  of  a 
straight  breakwater  at  an  angle  to  the  shoreline.  For  the  case  of  a  straight 
breakwater,  the  phase  function  S  lor  the  incident  wave,  the  phase  function 

S  for  the  reflected  waves  from  the  breakwater,  and  the  phase  (unction  H 

s  s 

(or  the  radiated  waves  generated  by  the  oscillatory  point  source  at  the  tip  of 
the  breakwater  can  he  expressed  respectively  as: 
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k  cos  0  <lx  +  K  v 

o' 


>  =  -  /  k  cos  8  dx  -  K  v 

s  /  o- 


=  -  J  k  cos  8  dx  +  j<ty 


K  -  k  sin  0  =  k  sin  0 
o  o  o 


K  -  k  sin  0  =  k  sin  0 
t  t  t 


k^  and  0  represent  the  wave  number  and  the  initial  angle  of  incidence  of  a 
radiated  wave  ray,  respectively.  The  numerical  model  developed  under  contract 
by  ini  based  on  the  uniformly  valid  asvmptot  it  theory  of  Liu,  Lozano,  and 
Pant  aza  ras  i  i ')  7‘i  )  was  mod  i  I  i  ed  to  conform  to  a  straight  breakwater  at  an  angle 
to  the  shoreline,  and  is  presented  in  Appendix  B  along  with  a  sample  output 
from  the  program.  This  numerii.il  model  served  as  the  basis  for  comparison 
with  t  tie  experiniiTit.il  results  of  the  present  study. 


I'.iVt1  s 

*  i  t  ti 

t  1  1  i  1  1  -If  *j-.v 

i  ’  f'l  •  h  :  1  .  J 

■  t  >  <  ■  ■  ■  < 

i  a 

add  l  t  , 

*ii  ,  tin 

i“  i  ‘ ' 

na  1,  t  e  .  - 

i  sled 

t  ha  t 

the  genes  at 

t  *  d  w.ivt-  hr 

1  gilt  s  si 

!•>(!'  ,  1  he 

1  a  I  ge 

enough 

■S  KJ 

l  1 1  a  t 

■  .11, a  1  1 

changes  m 

t  tie  tie  i  gli  l  s 

t  hr 

del  ei  la!' 

1  e  ,  a  rid 

at  Liu* 

same  t 

i  me 

l  tie 

tie  i  gli  t 

s  hull  1 

<1  tie 

sma  1  1  enough 

S'  t  1 1  •  1  t  t 

l  »!]i}'  tils-' 

ns  with 

i  i  lir  a  r 

t  tie.)  I  J 

e  s 

i  mild 

he  per- 

• 

!  (i  rim* 

.1.  1' 

re ! i mi na i A  i 

» •  S  l  <  ,  .Mid 

*  •  XlJ»  !  1  .lit 

• : .  1  a  1  re 

suits  o 

i  Hal e s 

t.  I 

dot)  /  , 

i  mli  - 

i,i  ted 

Lfiat 

tor  t!ie  are 

<1  of  lll.ljor 

i  rite  r'r  s 

1  ,  and 

lor  t  lie 

range 

of 

'..a 

oe  r  l  ,)d  s 

cons  i 

de  red 

pert  i  :ien t  , 

L  tie  .  t  v  1  f 

i  i  ex  ;  '«• , 

; men i a  1 

i*0 lid  i  L 

ions  shown 

in  7 

aide  1 

could  be  invest  ladled  with  height  changes  remaining  essentially  linear,  thus 
pe  mi  i t l  ing  irinp.iri.--mi:'  w  i  11,  t  iir< .  i  e  t  1 ,  i  i  development.,. 

■f  •  i  hi  r !  y i  x  w.ivr-he  i  ght  sms-irs  were  usi*<i  to  determine  the  wave 
heights  . i  1  mi g  t  tie  t  mi r  sei  t  i  mis  p.« r ■:  :  i  < *1  with  the  shore  1  i  lie  .uni  downw.i ve  of 


tin  *. , 

art 

i  c  a  i  l> 

I  f  J  K  Vs 

ate, 

f  t- 

i;;h:  >4 .  i  v.  i 

*.,  1  loll 

two  S 

ee  t  1 1  ms  ,  nine 

gages  a  1 ong  t w, 

St“t  t  It 

» 1 1  s 

,  .Mid 

U  :  1  1 

o  I  o  1  r 

lit 

o  gages 

ill  tile 

f.ifi  1 

egimi  lie. ir  tin 

‘  V 

ave  generator 

Idle  a* 

v’e  f 

age  ot 

l)t>  w 

JVr  h 

J  l 

X  h  t  s  i  •  • »  • 

,rded  , 

1 

ea,  h 

( >  t  the  two  in  t 

win 

gages  was 

<.  '  <ri  s  i  « 

ier 

ed  to 

he  t  h 

*■  :  up 

lit 

. ,  v  (  •  \  j , . 

:  ud  i  f  1  r 

'  Mil 

t  i  1  <* 

wave  genera t o r 

ihe  average 

of  <)(> 

.  1 

v  e  !  1  e  i 

gl  t 

r  .  j  i 

.  j. . 

1  ,1  lie 

f  -  Mil.*  1  II 

J  11 

m  >  -4 

gages  was  normal 

1  /.  d  to  t  il  l  s 

1  1 1  p  II  t 

■  1 

\*  e  1  a  •  i 

gill  . 

A  i  1 

*,n 

c’j  v.  t 

■if  rev 

« *  r 

■led  s 

i mu i t  a  neons  1 y 

V*  ) 

t  ti  ea  ch  Midi- 

V  1  ( 1  i  i .  1  . 

1  r 

e , .  d  1  1 1  g 

at  a 

1  i  d  J 

t  J 

slat  1  .  !  1 : 

•  (.  o  1 1  s  I 

s  1 

lllg  o 

1  Lite  average 

ot 

t  hose 

90  waves 

Ten 

repe 

a  t  a  b  i 

1  i 

! y  rep] i ; 

.  a t  i on 

le 

sts  were  conducted 

under  identical 

test  ; 

on 

d  i  t  i  mi 

S  to 

de  tin 

t  * 

t  he  Vail, 

1 1 )  i  1  i  t  v 

O 

f  t  tie 

me  a s u remeuts 

(e 

xpe  r imenta  1 

S  l  ,  1  t  t  * 

'  r  ) 

.  1  tie 

slat 

ist  i  •  ■ 

a  i 

Mi'  J  :  Ji  i  <  • 

of  tin 

V 

a  r  I  a  1  > 

l 1  i t  y  wa s  t  he 

sijuare  root  ot 

til"  V. 

1  ('  1 

JIM.  f  , 

or  t  ti 

r  :*,  t  i 

mf 

.Jr! 

it  1  ■  ■! , 

1  ' 

he  va 

r  1  a  lice  IS  deli 

Me* 

d  as  t  he  sum 

of  111- 

s 

jUJ  ifs 

*■1  i 

ti.  d. 

■ : 

■ 1  * 1  *  "■  ■ 

*•"  !'  .»• 

d  i 

V  1  dll.  i 

1  od'Se  1  V.I  t  loll 

1 1 

■  mu  tli.'ir  ave  i 

•  i  V,  e  ,  1 

I  :  v 

1  de  i  lr 

•  ■ 

i 

It  J  ! .  i  f !»  ’  t 

■■'*  •  o 

lllllt'e  !  "  i 

1  dec  |  a  t  es  .  (  111,- 

Standard  dev 

.  1  !  loll 

’a  ■  1 

S  1  i  S  1 J .  1 

lly  ■  • 

i.i*  . 

t 

j  .  <  1  .  . 

n*  ' 

1  1  ge  va  III.  O  t 

t  h 

e  observ.it  ion 

-4  -4 

idle 

W-lVf 

“  h  f  1  j.' 

i.i 

d  J  1  . 1  \%  1  *  i 

e  .  .J .  1  , 

1  ID 

f  d  VI 

til  the  wave  s  e 

1 1  s 

o  r  s  p<  is  i  L  i  ouei 

at  .'-I 

i 

inter*, 

t  i  s  s 

tail  I 

in: 

a  :  1  lie  1 

I  * 

i  .Hid 

extending  j 1 o 

n.u 

t  lie  sec  t  1  oils 

pa.  i  1  1 

e  i 

v.  1  l  h 

il"l  '  1 

:  n 

-  1  ■  W  1 1 W  • , '. 

f  '  *  I.I 

tie  IW 

t*  j  k'vN ,i  t  i*  r' .  !  o 

i"“' 

mV i dr  hr  t  t  r  r 

'lot  l  I.  1 

i  i  ■ 

Ol,  o| 

(lie  W 

■  1  V  t  •  ) 

*- 

!-■  ’!"■  - 

'*  * 

*d  1 

set  a  distance 

o 

1  ]  it  along 

tin*  - 

•  i  ! 

I"1*  P 

t  1  1  i  i 

-  ;  ,  . 

i  i  i 

.  :.*  .,*,  , 

•  :  ,  o 

1  : 

I  :  '  •  1  j 

V  f  \M  If  rep-si 

l  1  ■ 

■•led  at  J-  t  I 

* 

1  • 

?  f  .  f 

• . 1  iii' 

Is.. 

1  Mr 

d  1  li  t  o  i  ma  t  ion 

j  J 

i  -  t  ;  I  1 1 1  e  i  - 

i.  j  i  .  >  ; 

r 

i  :,*, 

!  JIM  !■ 

of  1 

I  *  If*  ’ll,  t 

» . , .  |ir. 

,!.* 

■'  *  1  t  t  I  ' 

along  two  sec 

l  1  < 

oils  and  I'M  a 

■  !  i  :  .  f.  -  1  .  : 

o  !  1  8 

;  t  i 

:  (  -;ii  ! 

hr' 

’  1  .  I  dr’..d‘ 

*•  ,  1  . 

:  i  r 

\  in- 

i'eili.1  1  u  i  n g  two 

SO 

i  t  Mills. 

1  ! 


B 


re 


« 


I  « 


,.1-ic  silt  1  i  ;  lent  I-.  mii  r  !  1  m  :  i;  that  the  u  :  v.  we  i  •  ,  •  '  :  .  .  ■„  i  t  -  -, 

M!icmi::l  i  eg  ion  .sin  mi  :n  i  :  1 1  r  •  ■  !.  i  lie  degree  I  I  .  !i<m  i  i  t  v  i  .  :  ;  .  •  \ 
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t  lie  average,  97.  Sd  pen  ent  el  the  ent-rpv  n|  tin-  watt  t  -  ■  i  r  *  -  x  i  s.  t  at  t  f  i 
eiated  pe  r  i  od  i  t  mnlamei, !  a  I  !  r-  ,  ,ei)  t  -  ‘ft  tie  inu  in  up  .  •  v  p<  t  . 

energy,  d  .  <9  pel  ,  rut  ex  I  si  s  at  tr  !  i  i  li  i  i  m-  a.  :  .  '  ’:e  .  :  a  t  i  :  1  -  i 

wave  energy  throughout  the  i;  •  n  !  !ur  tlie.se  represent  it  ive  t  ■  pei  une.'it 
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are  detectable  as  the  wav*  propagates  shoreward  past  the  a  i  i  •  ■  :  i,  ■  : 
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I  ri"  piency  to  higher  ha  mum  n.-,  i  nd  along  the  sei  :  ■  >>n  at  a  'titan-  . .  t 
t  com  the  shoreline  only  hd.uS  p.  rcent  exists  t  mi  i  i  e  av  age  i  a’  I  he 
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Figure  14.  Effects  of  nonlinearities  (initial  wave  height  and 
shoaling  effects)  on  percent  of  total  energy  in  fundamental  pe¬ 
riod,  wave  period  =  0.75  sec 


Figure  15.  Effects  of  non  1 i nea r i t ies  (initial  wave  height  and 
shoaling  effects)  on  percent  of  total  energy  in  fundamental  pe 
riod,  wave  period  =  1.00  sec 


K  i  gut  >•  !  t>  Kftcits  « .  t  n<  n !  1  ikm  r  1 1  1  es  (initial  wave  height  and 
sti".-1  :ry  ritri  isl  "ii  percent  of  total  energy  in  fundamental  pe- 
tied,  wave  period  =  !  .  50  sec 


w a ve -  He i gh  t  Da t  a 

*  7  U  iv  •  patterns  •  •  :  tug  t  lie  conditions  tested  in  this  experimental 
■  i  -  *•  ■. (roi.ii  in  I’ti.a  I-’t  wave-height  amplifnation  coefficients,  H/H 

I"!  all  i  ond  i  t  i  oiis  testei  a  r  •>  shown  i  u  Plates  l-.lt'.  Each  of  these  plates  con¬ 
stitutes  a  iei  t  ion  pai  t I  lei  with  the  shoreline  at  distances  from  the  shoreline 
>t  h ,  H,  10,  or  Id  tt  I h e  corresponding  si atter  of  the  experimental  data  is 
indiiated  hy  t  lie  error  Pars  of  p  !  us -o  r-in  i  nils  one  standard  deviation.  Also 
-hewn  in  these  p|  j  t  e  s  )r,.  i he  wave-height  amplification  due  to  diffraction 
i  lone  tas  i I  the  Ims  pi  were  horizontal  beyond  the  structure)  and  the  wave- 
:>•  i  ght  amp  I  i  t  i .  it  a ;  ini  1 1  at  e  d  hv  the  uniformly  v.,  1  i  d  asymptotic  theory. 

Study  Results 

48 .  for  small  wave  periods,  small  initial  wave  heights,  and  deep  within 
he  shadow  zone  (sections  near  the  shoreline),  both  the  diffraction  theory  and 
:  he  uniformly  valid  asymptotic  theory  predict  wave-height  amplifications  that 


compare  consistently  welL  with  the  experimental  data  Nearer  the  tip  of  the 
breakwater  the  diffraction  theory  tends  to  diverge  rapidly  from  the  experi¬ 
mental  results,  and  the  uniformly  valid  asymptotic  theory  more  nearly  approxi¬ 
mates  the  experiments.  As  wave  period  increases,  the  deviation  of  the  diffrac¬ 
tion  theory  becomes  more  pronounced  at  all  locations  except  very  near  the 
breakwater;  however,  the  uniformly  valid  asymptotic  theory  continues  to 
closely  approximate  the  data  from  the  experiments  in  all  regions.  The  greatest 
deviation  of  the  uniformly  valid  asymptotic  theory  from  the  experiments  occurs 
outside  the  breakwater  shadow  zone  in  the  area  of  the  asymptotic  undulations 
of  the  wave-height  amplification  factor.  This  can  be  attributed  to  reflection 
of  longer  period  wave  energy  from  the  experimental  beach  breaker  zone  not  ac¬ 
counted  for  in  the  analytical  development,  since  the  uniformly  valid  asymp¬ 
totic  theory  consistently  underpredicts  the  wave-height  amplification  in  this 
area  . 

49.  In  general,  it  can  be  concluded  that  the  uniformly  valid  asymptotic 
theory  is  superior  to  diffraction  theory  alone  for  estimating  wave  heights 
downcoast  of  nearshore  structures  subjected  to  combined  refraction  and  diffrac¬ 
tion.  Particularly  for  longer  period  waves  and  for  the  region  near  the  struc¬ 
tures  where  scour  and  erosion  are  known  to  frequently  occur,  this  theory  pro¬ 
vides  an  estimation  that  consistently  approximates  the  results  of  this  experi¬ 
mental  study  and  that  is  significantly  better  than  diffraction  theory.  At  the 
same  time,  it  appears  desirable  to  incorporate  into  this  theory  a  degree  of 
nonlinearity  that  is  not  presently  available. 

50.  Because  the  uniformly  valid  asymptotic  theory  compares  favorably 
with  the  results  of  this  experimental  study  (four  sections  parallel  with  the 
shoreline  for  three  wave  periods  with  three  incident  wave  heights  each) ,  the 
numerical  program  developed  from  this  theory  and  modified  to  fit  the  straight 
breakwater  geometry  was  applied  to  the  entire  shadow  zone  region  of  the  struc¬ 
ture.  These  results  are  presented  in  Plates  57-50  for  distances  parallel  to 
shoreline  of  1  ft  through  ! h  ft,  respectively,  in  increments  of  1-ft  displace¬ 
ment  seaward  toward  the  tip  of  the  breakwater  (which  was  located  15  ft  from 
the  shoreline).  The  three  wave  periods  utilized  in  the  experimental  phase  of 
this  study  (0.75,  1.00,  and  1.50  sec)  are  shown  in  these  plates  as  the  inde¬ 
pendent  parameters.  It  is  apparent  near  the  shoreline  that  the  longer  period 
waves  allow  for  the  greater  wave-height  amplification  at  all  locations  (par¬ 
ticularly  within  the  shadow  zone).  With  increase  in  distance  from  the 
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short-line,  the  longer  period  waves  continue  to  exhibit  a  greater  wave-height 
amplification  in  the  shadow  /one;  however,  outside  the  shadow  zone,  the  peri 
effect  becomes  less  apparent  and  eventually  becomes  entirely  obscured  as  the 
undulations  of  the  different  periods  engulf  (overlap)  ea< h  other. 


PART  VI:  CONCLUSIONS 


51.  Based  on  results  from  the  three-dimensional  experimental  investiga¬ 
tion,  reported  herein,  comparisons  with  existing  diffraction  theory,  and  a  new 
uniformly  valid  asymptotic  theory  reported  herein,  it  is  concluded  that: 

a.  For  short  wave  periods,  small  initial  wave  heights  (linear 
waves),  and  locations  well  within  the  shadow  zone  of  the  struc¬ 
ture,  both  the  diffraction  theory  and  the  uniformly  valid  as¬ 
ymptotic  theory  predict  wave-height  amplifications  that  compare 
consistently  well  with  the  experimental  data.  However,  nearer 
the  tip  of  the  structure,  diffraction  theory  tends  to  diverge 
rapidly  from  the  experimental  results,  and  the  uniformly  valid 
asymptotic  theory  more  nearly  approximates  the  experiments. 

b.  As  the  wave  period  increases,  the  deviation  of  the  diffraction 
theory  becomes  more  pronounced  at  all  locations  except  very 
near  the  structure;  however,  the  uniformly  valid  asymptotic 
theory  continues  to  closely  approximate  the  data  from  the  ex¬ 
periments  in  ail  regions. 

c.  The  uniformly  valid  asymptotic  theory  is  superior  to  diffrac¬ 
tion  theory  for  estimating  wave  heights  downcoast  of  nearshore 
structures  subjected  to  combined  refraction  and  diffraction. 
Particularly  for  longer  period  waves  and  for  the  region  near 
the  structures  where  scour  and  erosion  are  known  to  occur,  this 
theory  consistently  approximates  the  results  of  this  experi¬ 
mental  study. 
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table  1 


Experimental  Condi t ions  Tested 
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Wave  Height 

s  ,  H  (it) 
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near  Wave  Generator 
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Period,  T 

(  sec  ) 

0.7') 

1 . 00  _ 

1  .  50 

0.106 

0.102 

0.070 

0.  168 

0.130 

0.106 

0.218 

0.211 

0  .  la’ 

Note*:  All  experimental  test  conditions  were  repli- 

cat  ed  10  t i mes . 


Distribution  of  Wave  Knergy  Throughout  the  Kxpe  r  1  merit  a  1  Facility,  iO- 
Incident  Direction,  Period  -  0.70  sec,  (Kean  Wave  Height  -  (J  .  1  (J  t>  it 
tiO-deg  breakwater  Angle  with  Shoreline 
Percent  of  Total  Knergv  and  Phase  Angle  (deg),  Arrangement  Nh  .  J 


'I  able 


Distr  ihut  1  oil  ot  Wav-*  Energy  Ihroughout  t  fi.  •  Exj'er  i  wont  a  !  Facility,  3<)-<ieg 
Incident  Hi  feci  ion,  i'eno.l  -  0 .  75  sec,  Ocean  Wave  Height  =  O.ibH  it 
fiO-det;  Breakwater  Angie  with  Shoi.*l  ini* 

Percent  <>i  i'eta  I  Energy  and  Phase  Angie  (deg).  Arrangement  S'- .  1 
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0.01  (+76) 

n 

99. 29 (  +  1  34) 

0 . b  7 (-1 20) 

0.02 (+48) 

0.02 (+113) 

0.00 (+33) 

14 

99. 33 (-142) 

0.64 (+104) 

0.01  (-9) 

0.01  (-144) 

0 . 00 (—173) 

I  15 
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0.26  (-63) 
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0.01 (-15) 
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Distribut  ion  of  Wave  Energy  Throughout  the  Experimental  futility,  ;0-deg 
Incident  Direction,  Period  -  0 .  75  sec.  Ocean  Wave  Height  -  0 .  J.  1 8  it 


Pe 

60- 

rcent  of  Total 

deg  H  re  i kwu  t  e r 

Energy  and  Pha 

Angle  with  Shoreline 

use  Angle  (deg),  Arrangeme 

nt  No,  ^ 

Fundamental 

First 

Second 

Third 

Fourth. 

Period 

Harmonic 

Harmonic 

Ha  rmon ic 

Ha  rnon ic 

Gage 

0.75  sec 

0.38  sec 

0.25  sec 

0.19  sec 

0.15  sec 

1* 

94.49(+93) 

4 . 86 (+147) 

0.61 (-138) 

0.03  (-51 ) 

0.02 (+24) 

2* 

98 . 95 (+143) 

0.84  (-7 ) 

0 . 20  (+1  76 ) 

0.01  1+60) 

O.'JOi'  +  l  1 3 ) 

3 

93 .60  (  +  162) 

l .26 (+75) 

0.12 (-149) 

0.02  (+47 ) 

0. nO  (  +  146) 

4 

99.37 (-1 72) 

0.53 (-29 ) 

0.09  (+123) 

0.01 (-30) 

0.00 '+100 ) 

5 

95.55 (-73) 

4.25 (-134) 

0.16 (-179) 

0.02 (-161 ) 

0.02  (~“9) 

6 

95.43(+38) 

4. 32 (+76) 

0.1 8 (+117) 

0.01  (+85) 

0.01  (  +  33) 

7 

96 . 69 (+1 58) 

3.12 (-54) 

0.19  (  +  127) 

0 . 00 (-1 7 ) 

0.00  (  +  39) 

3 

97. 19 (-98) 

2.76 (+1 75) 

0.05  (  +  1 18) 

0.00 (-1 20) 

0.01  (-134) 

9 

96 . 34 ( +20) 

3.35 (+41 ) 

0.28  (+56) 

0.02  (  +  78) 

0.01  (  +  53) 

10 

96 . SO (  +  1 3 1 ) 

3.32 (-95) 

0.15 (+57 ) 

0 . 00  (-1 71) 

0.00  (+1 74) 

11 

93.13(147) 

1.83 (-1 38) 

0.04  (-17) 

0.01  (+64) 

0.01  (  +  179) 

12 

99.52 (+35) 

0.32 (+13) 

0.14  (+86) 

0.01  (  +  36) 

0.01  (  +  76) 

13 

99.09(+97) 

0 . 20 (+1 55) 

0 .09 (+108) 

0.01 (+77) 

0.01 (+76) 

14 

96 . 33  (  +  170) 

3.43  (-25) 

0. 24 (+139) 

0.01  (-5) 

0.00 (+123) 

15 

94.48  (-93) 

5.31 (-170' 

0. 19(+120) 

0.00 (+13) 

0.01  (-137) 

16 

92 .46  (  +  38) 

6.81 (+69) 

0.63  (+91 ) 

0.07  (+107) 

0.03  (  +  77) 

17 

94. 57 (+146) 

5.12  (-75) 

0.31  (+73) 

0.00  (-1 18) 

0.00  (-42) 

13 

95.25  (-102 ) 

4 . 59 ( +  1 65 ) 

0.15  (+70) 

0.01  (-60) 

0.00 (-14 7) 

19 

97.71  (  +  76) 

2.13  (+95) 

0.11  (  +  39) 

0 . 04 ( +93 ) 

0.01  (+92) 

20 

99 . 99 (+1 7 ) 

0.01 (-1 40) 

0.00(-120) 

0.00(+46) 

0.00  (+26) 

21 

99.62 (+21 ) 

0.36 (+35) 

0.00(-10) 

0.01 (+40) 

0.00  (  +  38) 

22 

95. 75  (+81  ) 

4.05 (+1 37) 

0.15  (  +  157) 

0.03(+77) 

0.01  (  +  79) 

23 

92 . 23 (+1 56> 

7.23 (-56) 

0. 54(+100) 

0.01  (-91) 

0.00 (+43) 

24 

90.02  (-92) 

9. 59 (+1 76) 

0.36  (  +  84) 

0.02  (-56) 

0.01 (-129) 

25 

84.85  (+2  7 ) 

12,h5(+46) 

2.06 (+64) 

0.33  (+76) 

0. 1 0 (+74 ) 

26 

89 . 60 (+1 53 ) 

9 . 3  5 ( —6 1 ) 

0.99 (+94) 

0.06  (-129) 

0.01  (  +  7) 

27 

9 1 . 80 (-103) 

7 . 94 (+160) 

0.23(+60) 

0.03  (-64) 

0.01 (-142) 

28 

98.06 (+125) 

1.92 (-103) 

0.02  (+39) 

0.00  (  +  13) 

0.00 (+39) 

29 

98,77 (+35) 

1.18 (+0.9) 

0.03  (-23) 

0.01  (-33) 

0.01 (+3) 

30 

99.6  7  (-2  7) 

0.37 (-31) 

0.01  (+54) 

0.Q1 (-85) 

0.00 (-62) 

31 

98.77  (-12) 

1.19 (-31 ) 

0.04 (+30) 

0.01  (-11) 

0.00  (-4) 

32 

91.39 (+25) 

7. 80 (+4 2) 

0.70 (+55) 

0.10(+69) 

0.02 (+60) 

33 

8  7.77  (+82 ) 

1 1 . 39 (+16) 

0. 73  (-1 35) 

0.03  (-8) 

0.03 (+45) 

34 

83 . 33 (+1 75) 

14.65 (-20) 

1  .  78 (-*-143) 

0.22 (-43) 

0.03  (+103) 

35 

31 .49  (-76) 

1 6 . 60 (-1 64 ) 

1 . 75 (+1 11) 

0. 12  (-10) 

0.04  (-122) 

36 

80 . 89 (+52 ) 

15. 36 (+86) 

3.01 (+108) 

0. 59  (  +  1 14) 

0 .  1  5  (  + 1  !  7  ) 

Ocean  gage  ■ 
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Distribution  of  Wave  Energy  Throughout  the  Expor  i  merit  .1 1  Facility,  30-deg 
Incident  Direction,  Period  =  1.00  sec.  Ocean  Wave  Height  -  0.102  It 
60-deg  Breakwater  Angle  with  Shoreline 
Percent  of  Total  Energy  and  Phase  Angle  (deg),  Arrangement  No.  1 


Gage 

Fundamental 

Period 

1.00  sec 

Fi  rs  t 

:1a  rnonic 
0.50  sec 

Second 

Ha rmonic 
0.33  sec 

Thirl 

Ha  rmonic 
0.25  sec 

Fourth 

Hi  rmon ic 
0.20  sec 

1* 

99.  19 (  +  86) 

0 .  o  5  ( + 1 1  6  ) 

0.1 6 (+174) 

0.01  (-133) 

0.00  (-25) 

2  * 

99.44 (+122) 

0 . 5  3  ( -96 ) 

0.00 (-125 ) 

0.02 (+1 3 ) 

0.01  1  +  119) 

3 

99. b 4 (+163) 

0.  32  (  +  1  56) 

0 . 0  4  (  -  5  6  ) 

0.00 (+128) 

0.00  (+2 7) 

4 

99. on (  + 1  17) 

0.3  31  +  147) 

0 . 0  5  ( - 1 4 1 ) 

0.02 (+10) 

0.01  (  +  113) 

3 

93.07 (+176) 

1.91 (-13) 

0.01 (+146) 

0.00 (-147) 

0.00 (+51) 

6 

98. 40 (-103) 

1  .  57  (  +  165) 

0.01 (+74) 

0.00 (-116) 

0.00 (-1 30) 

7 

93.  o.3  (-2  5) 

1  . '  ’  7  (  -  5 .3 ) 

0.01 (-72) 

0 .00 (-171) 

0.00  (  +  38) 

3 

9 3. 2 3 (+5 7  ) 

1.67  (+93) 

0.03 (+1 1 1 ) 

0.01 (+101) 

0.00  (  +  1 16) 

9 

93.75(+125) 

1.231-1 19) 

0.01 (  +  54  ) 

0.00  (+68) 

0.00 (+106) 

10 

99. 16 (-156) 

0.91 (+32 > 

0.02 (-136) 

0.00  (+99) 

0.00(-l 36) 

11 

99. 13 (-15) 

0. 84  (-22 ) 

0.01 (-60) 

0.02 (+142) 

0.01  (  +  142) 

12 

93.1  3 (-2 7) 

1.34 (-50) 

0.03 (-98) 

0.00  (  +  105) 

0.00  (+65) 

13 

99. 50 (-6) 

0 . 4  9  ( + 1 9  ) 

0.00 (-169) 

0.01 (+179) 

0.00  (-142) 

14 

93.3  3  (+64) 

1.63  (  +  105) 

0.04 (+106) 

0.00  (+140) 

0.00(+79) 

15 

9  7 . 9  3  ( + 1 4  4  ) 

1 .95  (-84) 

0.07 (+34) 

O.UO  (  +  1 56) 

0.00 (+61 ) 

16 

97.41 (-130) 

2.51  (+94) 

0.0  9 (-49 ; 

0.00(-l 50) 

0.00(-144) 

17 

93.44 (-63) 

1.841-121) 

0.01 (-1 59) 

0.00 (-6 7) 

0.00 (-77) 

13 

93. 26 (+13) 

i  .  7"'  t  +30) 

0.03  1  +  35  ) 

0.00  (-0.5) 

0.00 (+22 ) 

19 

99.76 (-11 3) 

0 .  !  6  !  -  .  1  ) 

0  .  U  7  (  + 1  7  7  ) 

0.01  (  +  151 ) 

0.00 (-32) 

20 

99.93 (-1  71 ) 

1  .o.c+126) 

0.04 (+16) 

0.01  (-131.) 

0.00 (+94) 

21 

99.9] (+170) 

0.03 (+3) 

0.02 (-18) 

0.01(4146) 

0.00 (-15) 

22 

93.95 (-166) 

1  .  n  3  ( 1 1  ) 

0.01 (+54) 

0.0  u -106) 

0.00 (+73) 

23 

97.021-109) 

2.“3i+l 33) 

0 . 0  5  ( +  7 ) 

0.01  (-126) 

0 . 00 ( + 1 2  7 ) 

24 

96.261-40) 

3.6  3  1-38) 

0.  1 1 (-113) 

0.00  (-109) 

0.00 (-141) 

2  5 

92.33 (+46) 

6.61  1  +  79) 

0.50 (+118) 

0.04  (+133) 

0.00  (+152) 

2h 

93.6  ! if  1 2) ) 

9 . 99 ( — 1 23) 

0.  3  5  1-6  ) 

0.04 (+1 1 3) 

0.00 (-9?) 

27 

{  -  !  S  ) 

’  .  +  3  6.  ) 

0 . 1 6  (- 1  2  5 ) 

0.01  (+83) 

n.OO  (-9  1  ) 

28 

93.97  *93  1 

0 . n ,  (+162  ) 

0.06  (  +  140 ) 

0.01  ('-165) 

0.00 (-17) 

29 

9  7.30  ( +6  1  t 

‘.121+38) 

0.04  (  +  13 

0.03  (+62) 

0.01  ( +98 ) 

30 

9  7 . 2  ’  ( +  .  2  ) 

2.676+47) 

0.0 2  1  + 8  5 ) 

0 . 04  ( -6  ) 

0.01 (+24) 

31 

95.61 (+49) 

4.  4 '10  3) 

0 . 1 1  ( + 1 2  3) 

0.02 (+100) 

0.02 (+102) 

32 

94 . 3  3  (  +  1 05  ) 

4.51  1-165) 

0.20 (-70) 

0.01  (-6) 

0.00  (  +  59) 

33 

90 .94 (  +  1 72  ) 

3.231-36) 

0 . 6  '3  ( + 1  20.) 

0.03 (-7 7) 

0.02 (+103) 

34 

85.03 (-1 10) 

1 2 .  34  (  +  1  24  ) 

1  .87 (  +  2  ) 

0.24 (-1 12) 

0.03 (+1 54) 

35 

83 . 24 (-3 1 ) 

1  3. 7 3 (-82) 

2.5 3 (—121) 

0.44 (-153) 

0.06 (-1 6  8 ) 

3o 

33. 32 (+53) 

16. 19 (+90) 

1  .  34  (+1  2  3  ) 

0.1  5 (  +  159) 

0. OOf-164) 

Ore 

an  gage. 

- P- 
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Distribution  of  Wave  Energy  Throughout  the  Experimental  Facility,  50-deg 
Incident  Direction,  Period  =  1.00  sec,  Ocean  Wave  Height  =  0.102  ft 


f > ( i  —  d e Breakwater  Angle  with  Shoreline 


1 

ercenl  ot  Total 

Energy  and  PI 

ase  Angle  (deg) 

,  Arrangement 

\n.  2 

Fundamental 

First 

Second 

Third 

Fourth 

Period 

Harmonic 

Harmonic 

Ha  rmonic 

Ha  rnon i c 

Cage 

1.00  sec 

0.30  sec 

0.33  sec 

0.23  sec 

0.20  sec 

I  * 

0d  .  1  0  7 +J0  .) 

0. 73(+L12) 

0.13  (-141) 

0.01 (-81) 

0.00 (+1 23) 

2* 

dd. 33  (  +  1 13) 

0.  34 (-? 7 ) 

0.037-29) 

0.0  3  (-33) 

0.01  (  +  76) 

3 

9d.2 3(+147) 

0.73 (+114) 

0.03  (-140) 

0.01 (+111) 

0.00  7-7>l  ) 

A 

94.28 (+159 ) 

0.71  f-1  01  .) 

0.01  7+2  3.) 

0.00 (-1 73) 

0.00/-! ! 7) 

") 

OS. 21 (-1 38) 

1.731+31) 

0.04 (-38) 

0.00  (+36) 

0.00 (-44  1 

h 

o 2 . 4 17-53) 

1  .  3  3  (-94 ) 

0.03 (-120) 

0.00  (+156) 

0.07)  (-1  7  '  ) 

d  3 .0d (  +  29 ) 

l .  ,3  3  ( +  3  3 ) 

0 .03 (+37) 

0.00 (+71) 

0.00(+l  34  ) 

44.0  ’  7  +  104 ) 

1 . 9d  (-1 43 ) 

0.0!  (-41) 

0. 007+93 ) 

o.no 7  +  1  5  7  i 

'  e  i  . ' '  7  +  I  7  7 >  i 

0 . 9  4  (—  1  7  ) 

o.Ol (+1 72) 

0. 00  (-100  ) 

0.00  7+  <J  ) 

1  0 

'd.l’f-'d! 

1.11  (  +  1  36) 

0.017+32) 

1.007-1  In) 

o .no 7 -:  +  1 1 

1 1 

<h  .  3rf  f  -*  2  n ) 

1 . 4  7  ( +'3:3  ) 

n  .lot-id]) 

0.03  (-5  3) 

0.0  I  7_2s) 

1  2 

' >7.'K :  1-20.1 

2.04  7-  1  7 ) 

0.02  (  +  31 ) 

0 .  ()•  !  (  +  143) 

0.01 7  +  1 47  ) 

1  3 

dd .37 (-1  l  ) 

1 .  37 ( +4) 

0.02  (-134) 

0. 037-14] ) 

n.ni  (-pm 

14 

99 . 6  5  ( +  34  ) 

0.3 3 7+ 31) 

0 .  o;)  (+ 1 0  3 ) 

0 . 00 (-6 1 ) 

0.0"7-14)) 

13 

93.4i r+ii3) 

1.347-143) 

0 . ' 1 3 (-39) 

0.02  (  +  31) 

o.  io7  +  ;  7 ) 

If. 

d  h  .  2 1  •  7 - 1  3  3  i 

3 .  4  -  f  +  )  3  ) 

).  :  -i  )'■: 

0.01  (  +  47  ) 

r . oi  '  .  i 

1  r 

d 7. 20 f- 79) 

2.  7  4  '  -  i  71  ) 

o.."  el  )■,  ) 

O.IO  (-40  ) 

0. 

1  < 

9 7  .  1 M 7 -4  ) 

7  1 

i . : :  7  - 1 1  > 

0.0 1  7-d  4  ) 

'  )  .  -  » 

l'» 

9d.  3 ( )  ( —  !  3  3 ) 

0 . 4  { -  4  1  ) 

0.07 (+40 ) 

I) . ')  3  (-  i2  ) 

n . 1 1 !  (  -  :  " 

y;  i 

44.3  7  ( -  1 hd ) 

('..0  1  (+174) 

'.(",+73) 

1.0 3 (-1 1 3) 

1 ) .  (>  1  (  +  );  i 

21 

44.  fe,  f-1  7.'  > 

' .  2  4  ( —  4  l  > 

0.0]  7+]24) 

0.0  H  +  l  T  l 

0.01  7  +  '  1  ' 

1  I 

4  7 . 4  3  <  - !  7  '■  ' 

.  4  ■  ♦ /  -v  H  0 ) 

7i.o ;  i-4  3 1 

O.no (  +  4  4 ) 

n.ni  i  (_i  ns ) 

2  3 

4b  .(H  7-34  ) 

3.3)7-123) 

o.  1  1  (-1  7  7  ) 

0.00  (  +  104 l 

0.00  (-14  1 

24 

4  3  .  6  .  +  1  ^  ) 

3.4o 7  +2  3 ) 

<>.  3n  7  + 3b  ) 

0.0)7+43) 

0.0  1  (  +  OI  l 

2  5 

4  ] .  n  3  ( + 1  o ; 

b. 42 7-172) 

0.40 (-73) 

0 . 0  4  (  +7  >4  ) 

0.00(  +  l 4  j  ) 

2r> 

44.3 < : - ! ’ 1 1 

4,237-3) 

0  .  i  4  7  -  1  7  6  ) 

0.0  3  7 - S ) 

0.01  f-1  *4  ) 

?7 

'■>  s  .  (-'i  :  ) 
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Distr ihut ; on  of  Wave  Energy  Throughout  the  Exper imenta 1  Facility^  30-Teg 
Incident  Direction,  Period  -  1.00  see.  Ocean  Wave  He ight  =  0.130  ft 
t>0-deg  Breakwater  Angle  with  Shoreline 
Percent  of  Total  Energy  and  Phase  Angle  (deg) ,  Arrangement  No.  2 
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Gage 
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D) st r i hut i  cm  n I  Wjvc  Energy  Throughout  the  Experimental  Facility,  JO-deg 
Incident  Direction.  Period  r-  1.50  so',  Ocean  Wave  Height  -  0.070  It 
Mi-  leg  In- ikw.it'  i  Angle  with  Shoreline 
Percent  *  *  t  Int.ii  Energy  and  Phase-  Angle  (deg).  Arrangement  No.  2 


Gage 

Fundament  a  1 

Pe  r  1  n<! 
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Vertical  direction  in  Cartesian  coordinate  system,  ft 
Arbitrary  parameter,  dimensionless 

2 

Battles  wave  number  correction  factor  (a  +  a  )/k  a, 

. .  .  ,  xx  yy 

d imensioniess 

Local  water-surface  elevation  and  complex  wave  amplitude,  ft 
Arbitrary  wave  angle,  deg 
Incident  wave  angle,  deg 

Reflected  wave  angle  along  breakwater,  deg 

Radiated  wave  angle  from  oscillatory  point  source  at  breakwater  tip, 
deg 

1.  1  U  1  302034  ,  dimensionless 


A3 


I/Vt"  inlt  ive  Ski- tell 


r'CrcCCCCCCCCCCOCCCCCCCCCCCCCCCCeCCCCCCCCCCCCCCCCC 
“  r  l  C  U  l  A  T  F  S  THE  WAVE  AMPLITUDE/  AND  DIRECTION 
-•A'.MI-'N  AT  A  CROSS  SECTION  X  =  Xl.  THE  TIP 

«  W  *.  T  f  0  IS  LOCATED  AT  X  r  XO/  T  =  TO. 

I  -  VALID  A  S  r  K  A-  T  0  T  !  C  THEORY  USED  IN  THE  ANALYTICAL 
ASSUMES  ’‘.AT  The  ORE  AKWATER  COINCIDES  WITH  A 
*  !  '■  T  F  0  fROM  the  TIP  OF  the  BREAKWATER.  HOWEVER/ 
i  version  or  the  program  has  been  modified  to 
The  more  usual  plan-form  oe  shore-connected 

AT  AN  angle  T-s  THE  SHORELINE/  MM  A  STRAIGHT  LINE 
I  T,.r  ray  a;  the  1|P  is  DESIGNATED  AS  ABR. 

CCCCCCCCCf.  ccccccccccccccccccccccccccccccccccccccccc 

f '• ) / F I (?,?>, GAUSS' ?,R) 
r  -  •  ~  ,  »c  E  f  (  11  0  1  ,  /  -1U  (  1  1  r,  ) 

.  c  ’  i '  -  /  s  t>  e  ( 1 1  i > ,  y  r  e  <  1 1  m 


•  •  -  r  < :  -  x  j  >  /  > . 

,  Hf.(-)l)/?, 


. :  c :  :c  cccccclc  ccccc  ccccccccccccccccccc 

i  v  T  .  a  I  A  F  3l  l  C  -  S  -- - - 


A  !  a  1  D  . 


c  c  c  r  -s  c  c  c  c "  c  '  t*.  c  r  c  c  c  c  ' .  c  c  c  c  c  c  c  c  c  c  c  c  c  c  c  c 


6  6  '  4  *  T  \  *  r  \  C  .  ;  -  |.u  s  »0 «  F  C  0  N  N  E  C  T  E  D  BPMk 

i  j  a  !  r  »  •  ' 

4  f?  I  i  c  <  1  ) 

6  *  f  )  R  *  A  r  ;  '  i  '  /  '  ‘J  J  v  f  ')  I  (  A  t  l  -I  I  t  w-  A  T  I  c  N  0  v  16  POINT  GAUSS 

1  0  ft  ;  Sr  A  T  t  C  • 

i  *  i  r,  .  l  !  .  *  o .  :■  (  i  '  ,  '  i  > 

.'1  f  -*•  ;  «  A  T  i  s  t  *  •.  •  i  s  i  *.  r  ;  u  \  •  > 


•  (  f  *  ’at  ‘  '•  1  /  1  M  S  m  *;  r  r-'  >  ;  •;  '  f  ft 1  ;• 

P  I  =  5 .  1  <•  1  ->  v  6  S  5  S  9 
O  =  l*;roi:/$f  "  i 
«.  >  =  i  \  t  :  K?  ./  $  r  r- ) 

K  S  F  A  ^ 

S  3  <  <  0  )  O  . 

A  3  4  (  <  n  )  =  ft  3  R  *  c  I  /  ■  0  n  . 

D">  \\\  '  ft  • =  i  /  1 

ccccccccrcccccc:rrccc:c:rfc':r::.:‘cccc(ccccrfcccccccccccccccccccccccccc 

. . T.»:  J  f  V  T  C  A  ^  r  f  .<  f  I  N  C  I  0  6  N  T  4  A  \J  F  DIRECT  ION - 

"  •  ■  '  •  •  .  '  •  (  •  r  c  c  r  •  C  C  f  r  c  c  c  c  c  c  c  c  c  c  c  c  c  c  c  c  c  c  c 

. .  -  r  .  ••  \t  »  7  T  *  A  i  •  ■  -  I  ‘Jr  T‘*L  AMOUNT  OF  COMPUTED  OUTPJ  T-- 


•.  L  «'  CC  r  i  CCCCCCCCCCrCCCCCCCCCCCC 

ST  r  '  -  ■*  A  r  ♦  ’ ,  \  i  *.  •  ,  1 

-••Me  (1  #6 :  i 

6?  >  )  R  -  A  i  •  *  ■  *  ,  '  ft  ■>  » i  •  A  f  !  N  f  .  P  C  R  I  0  D  X  OF  TIP 

1  x  "■  t  !  i' 


a  •'  r  a  *:  f  R  F  QUt  N  C  r  - 


«  1  -2S-C** 

*5  4. 


N  .  -  -  ►  c  at  f  ap  HFLO  (DEEP  4  A  T  F  R  >  • 
N  J  v  3  E  R  AT  T  I  3  — - - 


■;  *  T  A  *jm  <  P  *  M  )  ♦  C  *  p  •  H  /  (  C  C  S  H  (  P  *  H  )  ) 


I  t  |  p  .  -i  ;  ) 


.Ll.  3.005)  GO  10  5 


/  .  m  I  >  *  1  5  0  .  /  p  I 

;  z  ->  )  -  A  T  •  r'  !  /  1  S  0  . 

i 

'  •  '  •  p  f  [Hr"  r  , - 

.  .  .  r  r  7 'it  fli'lt  T  ION  OF  HE  INCIDENT  WAVE,  SO 

r  ;  'MS  OF  ThE  REFLECTED  WAVES/  A»EF, - 

' ;  .  s  t  (  w  '  *  •  ?  -  r".  r  d  *  *  ^  >  *  <  -  i .  >  *  s  t  p 

.  -  ’  .  •  o  j  o  t  (  j  n  r  •  •  ?  -  s  c :  •  •  ■’  ) 


,31.  x  .  N  [  .  1  7  GO  TO  ?  0  0  1 

the  L  r  C  a  T  I  ,7  N  OF  THE  BREAKWATER  '  ) 

v-r(V7  9?['iAT:  r-COD0f>I'JATF,> 


N  >:<T  »  L 
;  ’  N\  - <  ~ - L 

■  .  ,  7  0  UI 

■  •;  r> .  .  -  f  '  i  /  (  f.  .  T  A  M*1  (  f  .  H  )  ♦  r-  •  >’  *  H  / 

■  ■■",/!>)  , _  T  .  3 .  IIS)  GO  10  555 


'  .  )  f.  7  1  A  A  c 

1  . .  ?  -  -  ;  '  .  .  -> )  .  .  '  .  7  T  L  •  s  c  0 


.  .  J  -  ,  ;..1>)*O.S‘STP 

«  ’  -  .  /  (  ■  -  *  )  *  v  0  W  <  N  f.  )  *  S  C  l 


;  ■  o1*«  )  .  f  ri  *  »  i’  )  »  ,  f  *  c  •  'r  T  3 

r  (  .  1  *  •  :  f  s  #  s  ’  r 


s  r  i  ,'!(  '  :7',s  ■"  t  i  h  I  Rf  rifiTE  0 
-'■■as  i ■  i  r  o  m  i; - 


{  i'  «  h  )  ♦  r,  •  p  *  h  / 


(  Iisn  (  o.H  )  ,,  2  ) 

If  !A0S((P1-P)/P)  .LT.  0.005)  GO  10  772 
P  =  P  1 

GO  TO  771 

SCD=P1*SIN(AREF(NN>) 

P  f  f  <  (  ')  '4 )  -  1  C  0 


(VI  .  ♦  o  . 

1  ) 

5 1 =SC  D 

(  V  N  .  r  Q  . 

<  ?.) 

S  A  2=  S  C  0 

(K  1  .  P  Q  . 

NN  > 

STPN=Xl-x3 

(  <  '  •  N  t  * 

NN) 

STPN=A9S((x1-x3)/(K1-NN)) 

<  \>  »l  c 

<  1  ) 

0C  =  1  . 

*NN  .51. 

K  1  ) 

DC  =  -1  . 

8  A  Y  I  I =  3.  5  •  s:  0 /S  OR  T  (  PI  *  •  2  -  S  C  D*  *?  >*  (-1  .  '  ■  C*ST»N 
PhA1Is'1.5*STpn*(-1.)*DC*SQRT(P1»*2-SC.**  1 


(  N  N 

.EQ.  K1> 

K  =  1 

(NN 

-NE  . 

<  1  ) 

X  =  (< 1 -NN ) 

(  < 

.LT.  3)  < 

=  -< 

(  .N  N 

•  EQ. 

K  1  ) 

X  =  X  1 

(  X 

.  EQ. 

X  X) 

GO  TO  888 

0  0  R  S  ?  L  =  1  *< 

If  (NN  .'IE.  <1>  X  =  ST  PN*L  »  DC*  X  J 
If  ( L  .EQ.  1 )  GO  TO  883 

R  A  Y  I  IsRAYlI-DC/SQRT(p1**2-SCD*»2)»0.5*STPN»SCD 
PHAII*PHA!I-DC»SQRT(P1**2-SCD»*2)*0.5*STPN 

H  *  I)  E  P  r  H  (  X  ) 

P1sP-<G*P*TANM(P*H)-f2)/(G*TANH(P*H)*G»P»H/ 
COSH (P*H) »*2) 


IF  (A8S(<P1-P)/P)  .LT.  3.005)  GO  TO  881 
»  =  P1 

GO  TO  83? 

IF  ((SCO/PI)**?  . ST.  1.)  GO  TO  887 
R  A  Y  I  I  =  R  A  Y I  I-0C/S9RT (PI *  «  2  -  S  C  D  **  2) *.5*STPN*SCD 
PHAII=PHAI1-0C*SQRT(P1.*2-SCD«*2)».5»STPN 
CONTINUE 

YRE(NN)=YBW(NN)*RAYIl 

SRE(NN)pS0(NN)*A3$(PhAII>*ABS(X3-X)/<X3-X)»SC0»< 

Y  RE  ( NN) - YBW ( N  N )  ) 

GO  TO  770 
YRE(NN)=YBU(NN) 

SRE(NN)*SO (NN) 

WRITE  (10*230)  NN.YRE (NN) *SRE (NN) 

CONT I NUE 
GO  TO  85 
00  099  J»1  ,NN 
Ns  NN - J* 1 
"XlJltl 

YRE  (  M  )  *  Y  R  E  ( M  M  1  )  +1  . 

SR  E  (  " >  *SR  E (MM1 ) *S  CD*  1  . 

RE  F< ( Nlssco 

WRITE  (  10*  200)  M.yRE  (M),SRE(1M> 

CONT I NUE 
S  1  =S  C  0 

- compute  wave  number  at  the  cross  section- 

H  SDE  P  TH  (X  1  ) 

P1*d-(G«,»TANH(P*H)-F2)/(G*TANH(P*H)*G*P»H/C0SH(P»H) 

*•?. ) 

IF  ( *35 ( ( P1-P) / P>  .LI.  0.005)  GO  TO  17 

P  a  P  1 

GO  TO  16 
WNCS«P1 

ACS*ARSIN(SC  I  / WNC  S ) 

- COMPUTE  AMPLITUDE  RATIO  DUE  TO  REFRACTION 

AMF*$QRT( A8$ (WNCS*COS (AI*PI /I 80.)/WNI/C0S(ARSIN(SCI /WN 
C5))/(1.*2.»WNCS*DEPTH(X1)/SINH(2.»WNCS»0EPTH(X1))))) 
AMF«AMF*SQRT(1.*?.*WNI/S|NH(2.*WNI>) 

WRITE  (10*65) 

FORMAT  (’O'. ’AMPLITUDE  RATIO  DUE  TO  REFRACTION  ONLY’) 
WRITE  (10*130)  AMF 
WRITE  (10. «A) 

FORMAT  (’O’.  ’ANGLE  A  OF  INCIDENT  RAY  AT  THE  TIP*) 


u.*?  i  T  C  (in, 130V  a  t 

C - COMPUTE  WAVE  MO'S  AT  GAUSS  QUAD.  X-POINT - 

00  0  <=1,1 

<  <  =  3 

XG  =  (X1*X2>/7.*<X1-X0>/?.*GAUSS<1,IC> 

7  H=o;  PTxOr,  ) 

R  p  1  =■  =  -  (  *  -  .  T  t.  \  ■"  (  P  •  X  )  -  f  ?  )  7  (  C-  *  T  A  MM  (P*M)*S*P*H/(C05H(P*H)) 

1  • ) 

IF  (  A  0  S  i  l  p  1  -  P  )  /  •"> )  .LI.  0.005)  GO  TO  9 

p  =  p  1 

G  3  10  5 

9  IF  ;<<-!)  in , 1 • , 5 

10  W  M  (  i  •  <  -  1)  =  >'■  1 

<  <  :=  y  V  +  1 

X  G  =  x  1  +  X  V-  X  C. 

G  n  TO  7 

11  AM ( 2  *  X  )  =? 1 
6  CONTINUE 

C - C  a  Lf.  :L  a  t  F  incident  w  A  v  £  R«Y  PHASE  INTEGRAL  FROM 

C - TIP  T3  x-SEC  TIOM - 

P  H  A  I  I  =  0  . 

DO  T  2  1=1/3 

PHAII=PM»I!»(C0TAM(ARSIM<SCI/WN(?*L-1)))+C0TAN(ARSIN( 

1  SCr/WM(7*l))))*SCI*GAuSS(?.L)*(x1-xO)/?. 

17  CONTINUE 

WRITE  ( 1 n , 61 ) 

61  FORINT  (  • 0  '  ,  •  ANGLE  a  (DEG)  Y  LOCATION  amp.  ratio  dir 
1  .OF  nave  S'  ) 

c - PREPARE  S  C  A  L  L  I  N  G  - - 

M  1  s  A  1  /  A  5 

M7  =A  ’  /A  s 

C - ENTER  SCANNING  LOOP - 


DO 

13  R  =  '  1 

,m2 

A  = 

AS*1' 

I  F 

(a  . :  o . 

or.) 

G  0 

T  0 

1  3 

I  F 

(A  .  t  Q  . 

1  m. 

V  GO 

TO 

1  3 

I  f 

(  A  .  E  3  . 

2  70. 

)  GO 

r  o 

1  5 

I  F 

(A  .50. 

ARP) 

G  0 

T  0 

1  3 

c . ---INITIALIZE 

°  -1  A  I  D  =  0  . 


R  A  Y  I  D  =  "i . 

S  C  D  =  A  N  T  •  5  I  N  (  A  .  p  I  /  1  s  o  .  ) 

ADCS-ARSIN(SCD/WNCS> 

C _ COMPUTE  DIFFRACTED  RAY  R  AND  Y  INTEGRALS  FROM 

C - TIP  TO  X  -S  EC  T  I  ON - 

DO  15  M=1/« 

p.|AID=PHAID»(COTAN<AR5IN(SCD/AN(2*N-1)))»COTAN(ARSIN< 

1  5  :D/*M  (?  *N>  >  >  >  *3CD«GAij5S  (  7.N  I*  (X  1-XO)  /  7. 

RAYID=PAYIDfCTAN(ARSIN(SC0ZWNC7*N-1)))HAN(ARSlNCSCD/ 

T  <TN(2*N))))*6AuSS(2/N)*(x1-xO)/2. 

15  CONTINUE 

Y  =  -1  .  •  R  AY  I  D 

C - CO"0  HE  - - 

R  =  A  R  5  (  “-IN  I  D )  ♦A0S(Y*5CD) 

f - COMPUTE  - - 

S  =  A  T  5  (PM A  I  I)  •A3S<X0-X  1)/  (XO-X  1  )  ♦  Y  *  S  C  I 


C 


607 


F  D  1  =  1  . 


C  ...  ;-u!E  5  rl 


I  F 

(  ft  3 

R  .  G 

T .  AT)  GO 

10  60  1 

I  t 

(  Y 

.  3  T  . 

Y  9  F  (  K  ?  )  ) 

5  0  =  S  R  E (X?  )  ♦  SK  2  * 

(  Y- Y  R  E (<  2  )  ) 

I  f 

(  v 

.  G  T  . 

r  RE  (  K2 ) ) 

R  E  <  =  5  <  2 

I  F 

(  Y 

.LT  . 

Y  R  f  (  1  )  ) 

S  0  =  5  R  f  (1)*51*(Y- 

YRF  ( 

1  )  ) 

I  f 

(  Y 

.lT  . 

f  R  F  (  T  >  > 

R  E  <  =  5  1 

I  F 

(  Y 

*  J  *  • 

ypf(  n  . 

AND.  Y  .LE.  YRE (X2) ) 

GO  TO  60? 

GO 

T  1 

6  0  3 

DO 

6  0^* 

i=  1 

I  F 

(  Y 

•  LF  . 

Y  R  E  (  I  )  ) 

GO  TO  605 

CONTINUE 

l  F 

(  Y 

.  *  0  . 

r  (  [  )  ) 

5  P  =  S  R  E  (  I  ) 

B(v 


604 

605 


631 


606 


lf  <  Y  .fa.  YRE  (  T  )  )  RF<=Rf  F<  (  J  , 

If  <Y  .LT.  Y  R  E (  I >  )  S9=(Y-YRE(  I)  )  .  (  S  R  E  ( I-1J-SRE  (J) ) 
/(YRECI-1  1-YRE  (I  )  )  *  S  RE (  I  ) 

If  (V  .LT.  YRE(I)>  R E <  =  ( y-yre  (I  )  >*(REF < ( 1-1 ) -REF< ( I > ) / 

<refk<i-1)-rffk<i))*rffk(]) 

GO  TO  631 

If  (V  .31.  YRE(D)  S3  =  SRE  (  1  )  f  SI  *(Y-YRE  (  1  )  ) 

If  (Y  .31.  Y  R  E  <  1  1  1  REK  =  S  1 

YRE(K?)>  S  8  =  S  R  E ( K  2  )  +  SK  2  * (  Y  -  Y  RE ( K  2  ) ) 

RF<=SK? 

-  6  rg  0  .  Y  .  L  E  .  Y  R  E  (  1  )  )  GO  TO  606 


I  f 
I  F 
I  f 
GO 

DO 

If 


607  CONTINUE 


.  3T  . 

.GT  . 

.L  T  . 

.LT. 

.GF  . 

633 

1*1, <2 

YRE(  I)  ) 


(  Y 
(  Y 
(  Y 
(> 

T  0 

607 

(Y  .SB. 


Y  R  E  (  <  7  )  ) 

Y  R  F  (  <  ?  ) 


GO  TO  608 


33 


34 

1 

2 
3 


638 

I  f 

(  Y 

.ED. 

YRE (  I )  )  S9=$RE {  I  > 

I  f 

(  Y 

•  E3  . 

Y  R  E (  I  )  )  REK  =  REFK(I> 

1 

If  (Y  .GT. 
/  (  YRE  (  I- 1 

YRE(I))  S0=(Y-YRE(I)).(SRE(I-1)-SRE(I>) 

> - YRE  ( I ) )  +  SRE  (  I  ) 

1 

If  (Y  .GT. 
/  (YR  E  < I -1 

YRE(D)  REA=(Y-YRE(I>)*(REFK(I-1)-REFKCI>> 
)-YRc(I))*REE<(I) 

- - Uh  I  Hr  1  fl/  THE  TAB  F  j  R  DIFFERENT  REGIONS 

603 

I  f 

(  A  T 

•  GT 

.  A 9 R )  GO  TO  633 

I  F 

(  A 

.GT  . 

AT)  E  0  1  =  -  1  . 

I  E 

(  A 

.GT. 

A  9  R  )  F  D  1  =  1  . 

I  E 

(  A 

.  GT  . 

A  9  R  )  F  0  2  *  1  . 

I  E 

(  A 

•  GT. 

AREF  (<0)  *  1  RO.  /PI)  FD  2*-  1  . 

GO 

T  0 

33 

633 

I  f 

(A  . 

GT. 

ARE F ( K 3 ) • ) 83./P I  .AND.  A  .LT.  A9R)  FD2*1. 

I  E 

(  A 

.GT  . 

A9R  .AND.  A  .LI.  AT)  F  D 1 = - 1  . 

- CALCULATE  THETA,  THETAB - 

TRETA=SORT (R-S) 

If  (  ( R-S? )  . L  T.  3. ) 
THET  A  9  =  S  0  R  T  (  R-S0J 

. . . COMPUTE 

ESS*SIN(S  > 

E$C=CO$CS) 

ES9$=S  I  N ( S  9 ) 

£  S9C  *  COS(  S  3) 
XX=THETA* • 2. 


GO  TO  13 

SINES  AND  COSINES  OF  S,  SB- 


FRESNEL  INTEGRALS  VIA  SU9R0UTINE  CS - 

DO  35  N  N  = 1  ,2 

'  •  u'- 

CALL  C  S  <C  ,  SS  ,  XX  ) 

F  I  (NN,  1)  =  C 

FI  ( N  »J  ,  2  >  =  S  S 

• 

X  X  = TH  E  T  A  3  *  *2 . 

35 

CONTI NUE 

AMPL  I TUD  E  RATI  0 - 

9876 


*FI(1,2))*fo1*(ES9C*fi<?,1)-£S0s*FI(?,2)>»fo2)»»2»((esc 

»ES9C*ESS+ES9S)/7.*(ESS*FI  (1  ,1  > ♦ E  S  C  •  F I  ( 1 ,2) ) • EDI  ♦ 
(ES3S*FI<2»1)yES3C*F[(2,?))»fd2)**2)/2.) 

AA 1  =  3 . 5*E0 1»  E I ( 1 , 1 ) 

AA2=0.5.ED1»FI(1,2) 

8 1 *3. 5*E02  *F 1(2,1) 

R?  =3. S*FD2*F I (2,2) 

C  1  *A  A  1  *  AA  1  *A  A  2»  AA  ? 

C  2*3 1 *9 1 *02*92 
Dl =AA1*91-AA2*32 
D  2  *9 1  *  A  A  2  -  A A  1*32 
C  S  S  =  F  S  C  •ESBC*ESS*ES3S 
SSS=£SS»ES0C-ESC*ES9S 

DSINYsC1*SCI-C2*R£S*(SCI-RFO*(D1*CSS-D2*SSS) 

If  <A8S(REK)  ,GT.  A3S (WNCS) )  GO  TO  9876 

DS INX=C1*S0RT  (WNCS* WNCS-SCI»SCI  > ♦ C 2 • SO R T ( WNC S • WNC S -R £ K * R£<  > 

*^(S3RT  (uNCS.WNCS-SCl.SC  I) -SORT  (WNCS*WNCS-RE>(*REK)).(D1*CSS-D2 

GO  TO  376S 

DS INX.S9RT (WNCS* WNCS-SCI* SCI). (Cl* Dl* CSS-02* SSS) 
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99.000 

27.062 

1.177 

190.216 

1  00.000 

26.225 

1.215 

189.477 

1 01 . 000 

25.428 

1.267 

190.000 

1 02.000 

24.667 

1.291 

189.693 

1 03.000 

23.940 

1  .  275 

188.969 

1 04 . 000 

23.245 

1 . 230 

188.713 

1 05.000 

22.578 

1  .  183 

188.959 

1  06 . 000 

21.938 

1 .160 

189.326 

1 07. 000 

21.323 

1.171 

189.507 

1 08.000 

20.731 

1  .  210 

189.524 

1 09. 000 

20.161 

1.257 

189.575 

1 10.000 

19.612 

1.295 

189.741 

111. 000 

19.082 

1  .  31  4 

189.802 

112. 000 

18.569 

1.310 

189.359 

1 1 3.000 

18.073 

1.285 

188.533 

1 1 4. 000 

17.594 

1.245 

188.436 

1 15.000 

17.129 

1.200 

189.547 

1 16.000 

16.678 

1.159 

189.947 

1 17.000 

1 6.240 

1.131 

189.135 

1 18.000 

15.815 

1.119 

190.225 

1 19.000 

1 5.402 

1.127 

190.862 

1 20.000 

1 5.000 

1.150 

189.415 

121. 000 

1 4.608 

1.186 

190.823 

1 22.000 

14.227 

1 .228 

189.951 

1 23. 000 

13.855 

1.271 

189.695 

1 24.000 

1  3 .492 

1  .  31  1 

190.316 

1 25.000 

13.138 

1  .  345 

188.971 

1 26.000 

12.792 

1  .  37  2 

190.147 

1 27.000 

12.454 

1  .  391 

188.600 

1 28.000 

12.124 

1  .400 

189.817 

1 29. 000 

11.800 

1.401 

188.351 

1 30.000 

1 1  .483 

1  .  395 

189.532 

131. 000 

1  1  .1 73 

1  .  380 

188. 207 

1 32.000 

10.869 

1  .  360 

189.276 

1 33.000 

10.571 

1.334 

1 33.285 

1 34 . 000 

10.278 

1  .  303 

188.920 

135. 000 

9.990 

1.269 

183. 700 

136. 000 

9.708 

1.231 

188.419 

1 37.000 

9.431 

1.192 

189.291 

1 38. 000 

9.158 

1.151 

138.141 

139.000 

8.889 

1.109 

139.427 

1 40. 000 

8.625 

1 . 066 

188. 742 

141.000 

8.365 

1  .024 

188. 708 

142. 000 

8.109 

0.982 

1 89.854 

143. 000 

7.856 

0.94  1 

138.394 

144.000 

7.607 

0.901 

189. 568 

1 45 . 000 

7 .362 

3.861 

139.911 

1 46.000 

7.119 

0.  824 

188.303 

1 47.000 

6.880 

0. 787 

190.308 

148.000 

6.644 

0.752 

139.963 

1 49.000 

6.4  10 

0 . 71  9 

1  R8 . 1  1  7 

1 50. 000 

6.180 

0.687 

190.849 
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1  51 . 000 

5.951 

0 . 65  7 

190.272 

is?. ono 

5.725 

0.628 

187.550 

153. 000 

5.502 

0.601 

190.958 

154.000 

5  .  281 

0.575 

191.175 

155. 000 

5.061 

0.55  1 

1 86.741 

156. 000 

4.844 

°  .  5  7  8 

189.970 

1  57.000 

4.629 

0.506 

1 92 . 775 

1 58.0QO 

4.416 

0.486 

1 86 . 959 

1 50. 0O0 

4.704 

0.466 

187.290 

i  6  0. ono 

3.994 

0. 44  8 

193.374 

i 6i . ono 

3.785 

0.43? 

190.259 

167.000 

3.578 

0.416 

184.335 

163.000 

3.372 

0.401 

19?.041 

1 64 . 000 

7.167 

0.  787 

195. 523 

165.000 

2.964 

0.374 

183.374 

1 66 . °33 

2.761 

0.781 

133.658 

167.  mo 

7.560 

3.350 

197.111 

1 63.000 

7.359 

3.  339 

196.803 

1  69.000 

2.160 

0.  32  9 

192. 096 

170.000 

1.961 

0.319 

195.426 

1  7 1 . 000 

1.763 

0.310 

1  99.444 

177.000 

1  .565 

0.  302 

198.447 

1 7  3 . non 

1.369 

0.294 

198. 702 

1  7  4.000 

1  .  1  72 

0.236 

193.672 

175. 000 

0.976 

0.279 

194.461 

1  7  6.000 

0.780 

n.  273 

195. 521 

1  77. 000 

0.885 

0.266 

201 . 1 90 

1 73. ODD 

0.3  90 

0.261 

195.506 

1 79.Q00 

0.195 

0.  255 

1 82. 298 

181. 000 

-  0 .195 

0.245 

201.580 

1 87.000 

-5  ,  U7 

0.241 

197.383 

183. non 

-0.585 

0.236 

186.856 

184.000 

-  0  .  7  8  : 

0  .  23? 

195.187 

1  8.5  .  030 

-0.076 

0.278 

205.808 

186.  mo 

-i  .1  77 

0.  ??  5 

205.254 

1 87.00° 

-1.369 

1.222 

205.018 

188.000 

-1.565 

0.  ’i  8 

206.284 

1  8  9  .  0 00 

-1.763 

0.215 

196.671 

i  )  n .  n  m 

-1.961 

0.213 

183. 365 

191. non 

-7.160 

0.210 

200.896 

1^7. 

-7.359 

3.20  7 

206. 106 

1  9  7  . 

-7.560 

0.205 

191 .716 

19  4.  I?:.' 

-7.761 

3.20  3 

138. 340 

195.  mo 

-7.964 

0.201 

205.698 

1  g  n  t  t 

-3.167 

0.199 

210. 720 

197. 0°0 

-7.37? 

0.197 

208.663 

193.000 

-7.578 

0.196 

712. 1 82 

1  9  9  .  o  °  f) 

-7.735 

0.194 

204.668 

? on. 030 

-3.994 

0.193 

130.818 

?  oi . ono 

-4.204 

0.19? 

189.190 

707.00 o 

-4.416 

0.191 

211.458 

703.000 

-4.629 

0.  190 

212.071 

204. 000 

-4.844 

0.189 

210.642 

205.000 

-5.061 

0.188 

2  1  3  .  7  9  2 

2 06.000 

-  5  .  ?  r,  1 

0  .  1  8 .8 

198.754 

2  07. 000 

-5.502 

0.18? 

1 86. 079 

2  OR. 00  0 

-5.725 

0.187 

208 . 287 

7  09. 000 

-5.951 

0.18? 

71 4 . 457 

Definition 

of  Program  Parameters 

Name 

Description 

'A 


A 

Scanning  angle 

“•  *•  V 
-.1 

ABR 

Angle  of  tip  of  breakwater  with  normal  to  shoreline 

ABW(llO) 

Local  angle  of  breakwater 

• 

ACS 

Angle  of  incidence  at  arbitrary  cross  section 

ADCS 

Angle  of  radiated  ray  at  arbitrary  cross  section 

AI 

Angle  of  incidence  at  infinity  (deep  water) 

AMF 

Amplitude  ratio  due  to  shoaling  effect 

• 

AMR 

Amplitude  ratio  due  to  combined  refraction  and  diffraction 

APA 

Local  angle  of  wave  propagation 

AREF (110) 

Angle  of  reflected  wave  along  the  breakwater 

'  V 

AS 

Scanning  angle  step 

• 

AT 

Angle  of  incidence  at  the  tip  of  the  breakwater 

Al 

Lower  limit  of  the  scanning  angle  (>90  deg) 

A2 

Upper  limit  of  the  scanning  angle  (<270  deg) 

CS 

Fresnel  integral  subroutine 

• 

DEPTH (X) 

Water  depth 

, 

ESBC 

Cos(SB)  (Cosine  of  phase  function  of  reflected  wave  ray) 

ESBS 

Sin(SB)  (Sine  of  phase  function  of  reflected  wave  ray) 

ESC 

Cos(S)  (Cosine  of  phase  function  of  incident  wave  ray) 

• 

ESS 

Sin(S)  (Sine  of  phase  function  of  incident  wave  ray) 

FI (2 ,2) 

Output  for  subroutine  CS 

G 

2 

Gravitational  acceleration  constant,  32.15184  ft/sec 

GAUSS (2, 8) 

Location  and  weight  factor  for  16-point  Gauss  Quadrature 

• 

H 

Water  depth 

)  ■  ■  . 

PHAID 

Integra]  part  of  phase  function  of  radiated  wave  ray 

•  -  *  f 

PHAII 

Integral  part  of  phase  function  of  incident  wave  ray 

' .  j 

PI 

3.  14159265359 

• 

B1 1 


Name 


RAY  ID 
RAY  1 1 


SO(llO) 
SRK( 1 10) 


THF.TA 
I'HKTAB 
WN( 16) 


YRVl 110) 
YRK ( 110) 


Phase  function  of  radiated  wave  ray 

Integral  part  of  the  location  of  a  radiated  wave  ray 
Integral  part  of  the  location  of  a  reflected  wave  ray 
Phase  function  of  incident  wave  ray 
Phase  function  of  reflected  wave  ray 

Snell's  constant  for  radiated  wave  ray  from  breakwater  tip 

Snell's  constant  tor  incident  wave  ray 

Phase  functions  of  incident  waves  along  the  breakwater 

Phase  function  of  reflected  wave  along  XI  cross  section 

Step  size  of  numerical  integration  along  a  ray  by  the  trapezoidal 
ru  1  e 

Wave  period  in  seconds 
Square  root  of  (R  -  S) 

Square  root  of  (R  -  SB) 

Wave  number  at  Gauss  Quadrature  points 
Wave  number  at  arbitrary  cross  section 
Wave  number  at  infinity  (deep  water) 

Wave  number  at  tip  of  breakwater 
X-coordinale  of  tip  of  breakwater 
X- coordinate  of  arbitrary  cross  section 
X-coordinate  of  end  of  tracing  region  (X2  >  X0 ) 

X-coordinate  of  Gauss  Quadrature  points 
Y-coord  i  mat  es  of  the-  breakwater 

Y-eoord i nut es  of  reflected  wave  rays  along  XI  arbitrary  cross 

S('(  t  1  oil 


•>  •  ,  w  , 
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